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ABSTRACT 


The  present  work  is  concerned  with  development  of  innovative  damage  mechanics  models 
for  constitutive  modeling  of  modem  composite  materials.  Special  attention  is  focused  on 
micromechanical  damage  theories  to  explain  and  model  behavior  of  composites  based  on 
micromechanics,  instead  of  using  phenomenological  models  with  arbitrary  thermodynamic 
potential  functions  and  arbitrary  damage  evolution  equations.  In  particular,  basic  studies  are 
performed  on  the  notion  of  isotropic  and  anisotropic  damage  variables  in  continuum  damage 
mechanics.  It  is  demonstrated  that  "isotropic  damage"  does  not  necessarily  imply  "scalar  dam¬ 
age"  representation  in  general.  Furthermore,  a  micromechanical  damage  model  is  presented  for 
uniaxially  reinforced  composites  with  interfacial  arc  microcracks.  All  microcracks  are  assumed 
to  occur  in  the  fiber/matrix  interfaces,  and  are  modeled  as  arc  microcracks.  Microcrack-induced 
strains  and  overall  compliances  are  analytically  derived  based  on  micromechanical  bi-material 
arc  crack  solutions.  Both  "stationary"  and  "process"  damage  models  are  given.  It  is  emphasized 
that  the  present  work  does  not  employ  any  arbitrary  (fitted)  "material  constants".  Finally,  a  novel 
three-dimensional  statistical  micromechanical  theory  is  proposed  to  investigate  the  nonlinear 
behavior  of  microcrack-weakened  brittle  solids.  The  macroscopic  stress-strain  relations  of  elas¬ 
tic  solids  with  interacting  microcracks  are  micromechanically  derived  by  taking  the  ensemble 
average  (over  all  possible  realizations)  and  volume  average.  The  proposed  statistical 
micromechanical  theory  accounts  for  random  microcrack  location,  distribution,  and  interaction 
without  using  Monte  Carlo  simulations.  The  theory  is  fundamentally  different  from  existing 
effective  medium  theories  and  deterministic  microcrack  interaction  method. 


RESEARCH  OBJECTIVES 


The  proposed  research  work  is  concerned  with  development  of  innovative  damage  mechan¬ 
ics  theories  and  constitutive  modeling  of  modem  composite  materials.  Micromechanically 
based  damage  theories  and  models  will  be  actively  pursued.  Particular  attention  will  be  focused 
on  the  construction  of  physical  underlying  damage  mechanisms  based  on  microstructural  con¬ 
siderations.  The  general  goal  is  to  physically  understand  and  mathematically  model  the  funda¬ 
mental  damage  mechanisms  at  the  local  level  so  that  we  can  simulate  randomly  oriented  distri¬ 
buted  microcracks  and  perform  failure  analysis  of  composite  solids  and  structures.  The  specific 
objectives  of  the  proposed  research  are  the  following: 

1.  Develop  proper  micromechanically  based  "damage  mode  tensors"  for  different  modes  of 
failure  mechanisms. 

2.  Develop  a  set  of  physical  damage  loading/unloading  fracture  criteria  for  microcracks  initia¬ 
tion,  growth,  coalescence  and  arrest. 

3.  Derive  micromechanically- based  evolution  equations  for  microcrack  growth. 

4.  Derive  progressively  damaged  material  properties  (complainces)  and  obtain  damaged 
anisotropic  multi-phase  stress-strain  constitutive  laws  based  on  micromechanics  solutions. 

5.  Develop  efficient  and  accurate  computational  constitutive  algorithms  for  the 
micromechanically-based  damage  models  proposed,  and  perform  experimental  verification 
of  the  models  obtained. 
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PART  I: 

On  the  Notion  of  Isotropic  and  Anisotropic  Damage 
Variables  in  Continuum  Damage  Mechanics 

LO.  Abstract 

The  present  work  analyzes  the  implication  and  limitation  of  some  “scalar  damage”  models.  In 
particular,  thermodynamic  potential  and  “effective  stress  concept”  are  re-examined.  It  is  demon¬ 
strated  that  “isotropic  damage”  does  not  necessarily  imply  “scalar  damage”  representation  in  gen¬ 
eral.  The  notion  of  isotropic  and  anisotropic  damage  variables  in  continuum  damage  mechanics  is 
then  discussed.  In  addition,  some  results  from  micromechanical  analyses  are  applied  to  show  the 
direct  relationship  between  the  fourth-order  damage  tensor  and  the  damage-induced  compliance 
tensor  characteristic  of  microcrack-weakened  brittle  matrix  composite  materials.  It  is  shown  that 
even  for  isotropic  damage  one  should  employ  an  isotropic  fourth-order  damage  tensor  (not  a  scalar 
damage  variable)  to  characterize  the  state  of  damage  in  materials,  in  accordance  with  the  effective 
stress  concept.  In  general,  however,  a  damage  tensor  is  anisotropic  and  should  be  derived  from 
micromechanical  damage  analysis  when  possible. 


X 

1.1.  Introduction 

Since  the  first  introduction  of  the  scalar  damage  concept  by  Kachanov  (1958)  and  Rabot- 
nov  (1963)  for  creep  of  metals,  continuum  damage  mechanics  has  become  an  emerging  field  of 
active  research.  Extensive  phenomenological  damage  models  were  proposed  in  the  literature  to 
describe  and  predict  macroscopic  constitutive  behavior  of  ductile  and  brittle  materials  containing 
distributed  microdefects  such  as  microvoids  and  microcracks.  We  refer  to  Lemaitre  (1984),  Kra- 
jcinovic  (1984,1986,1989),  Kachanov  (1986),  and  Ju  (1989a)  for  a  literature  review.  In  addition, 
micromechanical  “non-process”  (stationary)  and  “process”  (propagating)  damage  models  for  brit¬ 
tle  solids  with  many  microcracks  were  proposed  by,  e.g.,  Budiansky  and  O’Connel  (1976),  Horii 
and  Nemat-Nasser  (1983,1985),  Kachanov  (1987),  Krajcinovic  and  Fanella  (1986),  Sumarac  and 
Krajcinovic  (1987,1989),  Fanella  and  Krajcinovic  (1988),  Krajcinovic  and  Sumarac  (1989),  and  Ju 
(1989b). 

The  predictive  utility  of  a  damage  model  depends  heavily  on  its  particular  choice  of  a  “ damage 
variable ”  which  serves  as  a  macroscopic  approximation  in  describing  the  underlying  microme¬ 
chanical  processes  of  microdefects.  In  the  current  literature,  there  are  many  ways  to  phenomeno¬ 
logically  define  or  micromechanically  derive  damage  variables;  see,  Krajcinovic  (1989)  for  a  state- 
of-the-art  review.  In  particular,  scalar  damage  variables  were  widely  used  for  isotropic  or  one¬ 
dimensional  phenomenological  damage  models;  see,  e.g.,  Lemaitre  and  Chaboche  (1978,1985)  for 
excellent  work.  In  addition,  vectorial,  second-order  and  fourth-order  tensorial  damage  variables 
were  often  used  for  anisotropic  phenomenological  damage  models;  see,  e.g,  Krajcinovic  and  Fon- 
seka  (1981),  Vakulenko  and  Kachanov  (1971),  Kachanov  (1980),  Murakami  and  Ohno  (1981), 
Chaboche  (1979),  Cordebois  and  Sidoroff  (1979,1982),  and  Ju  (1989a)  for  interesting  examples. 
The  simplicity  of  a  scalar  damage  representation  is  indeed  very  attractive.  However,  a  scalar  dam¬ 
age  model  is  somewhat  of  limited  use  in  practice.  A  vectorial  damage  representation  is  appealing 
because  microcrack  areas  and  orientations  are  captured  and  a  vectorial  representation  is  a  direct 
extension  of  the  original  Kachanov’s  model.  Nevertheless,  a  vector  damage  model  to  some  ex¬ 
tent  suffers  microcrack  shape-indifference  and  tensor  operation  problems  related  to  general  stress 
transformation.  A  second-order  damage  tensor  representation  is  also  attractive  since  microcrack 
areas  and  orientations  are  captured,  too.  It  is,  however,  incapable  of  describing  general  anisotropy. 
Therefore,  as  noted  by  several  researchers,  an  appropriate  description  of  anisotropic  damage  gen¬ 
erally  involves  a  fourth-order  or  even  an  eighth-order  damage  tensor  representation.  In  practice. 
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an  eighth-order  damage  representation  is  perhaps  too  complicated.  A  fourth-order  damage  repre¬ 
sentation,  on  the  other  hand,  can  be  handled  reasonably  well  and  and  has  nice  correspondence  with 
the  fourth-order  overall  stiffness  tensor. 

In  spite  of  the  direct  intuition  that  “isotropic  damage”  necessarily  warrants  a  scalar  damage  rep¬ 
resentation  (see,  e.g.,  Lemaitre  and  Chaboche  (1978),  Lemaitre  (1984),  Mazars  (1984),  Kachanov 
(1986)),  this  is,  surprisingly,  not  always  true  from  the  viewpoint  of  the  effective  stress  concept  and 
thermodynamic  potential.  The  present  work  shows  that  isotropic  damage  actually  only  implies  that 
the  (fourth-order)  damage  tensor  is  isotropic ,  not  necessarily  reducing  to  a  scalar.  On  the  other 
hand,  anisotropic  damage  implies  that  damage  tensor  is  anisotropic.  Further,  the  implication  of 
scalar  damage  representation  upon  overall  elastic -damage  moduli  is  discussed  within  the  frame¬ 
work  of  the  effective  stress  concept.  Some  results  from  micromechanical  analyses  are  also  applied 
to  demonstrate  the  physical  meaning  of  the  fourth-order  damage  tensor  in  the  case  of  microcrack- 
weakened  brittle  solids. 
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1.2.  Implication  of  scalar  damage  variable 
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In  the  case  of  spatially  perfectly  randomly  distributed  microcracks  or  microvoids  in  all  direc¬ 
tions,  isotropic  damage  model  is  certainly  an  appropriate  choice.  In  particular,  in  a  scalar  damage 
model,  the  virgin  (undamaged)  material  is  often  assumed  to  be  linearly  elastic  and  isotropic  in  the 
absence  of  plastic  flow.  In  the  literature,  some  researchers  consider  the  following  rational  elastic- 
damage  potential  energy  based  on  the  physically  appealing  effective  stress  concept  (Lemaitre  and 
Chaboche  (1978),  Lemaitre  (1984),  Mazars  (1984),  etc.): 

J  (1  -  d)  ee :  C°  :  £e  (1) 

where  d  is  a  scalar  damage  variable  (ranging  from  0  to  1),  ce  is  the  elastic  strain  tensor,  C°  is 
the  virgin  linear  elasticity  tensor,  and  “  :  ”  signifies  tensor  contraction.  From  the  Clausius-Duhem 
inequality  for  isothermal  processes,  one  obtains  the  homogenized  stress  tr  : 

<r  =  (1  -  d)  C°  :  £e  (2) 

Hence,  it  is  clear  that  the  overall  elastic-damage  stiffness  tensor  C  takes  the  form: 

C  =  (l-d)C°  (3) 

such  that  <r  =  C  :  e'  .  As  a  consequence,  a  component  of  the  overall  elastic-damage  stiffness 
tensor  can  be  obtained  by  scaling  down  the  corresponding  component  of  the  virgin  elastic  stiffness 
tensor  by  a  common  factor  (1  —  d) .  Moreover,  the  effective  stress  &  is  often  defined  as  (within  the 
framework  of  isotropic  damage): 


Thus,  we  have  &  =  C° :  ee .  The  above  formulation  is  indeed  excellent  for  one-dimensional  prob¬ 
lems  and  is  entirely  consistent  with  Kachanov’s  (1958)  original  concept  of  damage  measure. 

In  a  three-dimensional  setting,  nevertheless,  Eqs.  (l)-(4)  are  less  attractive.  This  can  be  proved 
by  simply  noting  from  Eq.  (3)  that  damaged  bulk  and  shear  moduli  become  I<  =  (1  -  d)I<0  and 
G  =  (1  -  d)G0 ,  respectively,  with  K0  and  Ga  denoting  the  virgin  bulk  and  shear  moduli.  It  follows 
that  Gj K  =  G0/K0  .  Therefore,  the  above  scalar  damage  formulation  inherently  implies  that 
Poisson’s  ratio  always  remains  constant,  i.e.,  v  =  v„  .  This  is  somewhat  too  restrictive  and  not 
universal  even  among  isotropic  damage  processes.  This  anomaly  is  due  to  the  one-dimensional 
nature  of  Eq.  (4). 
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1.3.  Notion  of  isotropic  and  anisotropic  damage  variables 

As  discussed  in  the  previous  section,  “scalar  damage”  is  only  a  special  case  of  “isotropic  dam¬ 
age”.  “Isotropic  damage”  only  states  that  the  fourth-order  damage  tensor  D  is  isotropic.  In  addition, 
isotropic  damage  does  not  require  the  virgin  elasticity  tensor  C°  to  be  isotropic.  It  is  emphasized 
that  isotropic  damage  will  preserve  the  directional  characteristics  of  the  virgin  elasticity  tensor  C° . 
That  is,  if  C°  is  isotropic,  then  the  elastic-damage  stiffness  C  is  also  isotropic.  On  the  other  hand,  if 
C°  is  orthotropic  (or  anisotropic),  then  C  is  also  orthotropic  (or  anisotropic).  By  contrast,  general 
anisotropic  damage  will  change  the  directional  properties  of  C°. 

The  anomaly  in  the  aforementioned  scalar  damage  formulation  can  be  easily  resolved  by  re¬ 
examining  the  effective  stress  concept  and  the  elastic-damage  energy  potential.  In  particular,  the 
following  formulation  is  free  of  the  previous  anomaly  in  three-dimension: 

V  =  \  £e :  [M(D)  •  C°] :  ce  (5) 

where  M(D)  is  a  fourth-order  tensorial  function  of  the  fourth-order  damage  tensor  D .  The  Clausius- 
Duhem  inequality  then  renders 

tr  =  [M(D)  •  C°]  :  ee  (6) 

According  to  the  effective  stress  concept,  the  effective  stress  can  be  expressed  as  a  =  C° :  ee  . 
Consequently,  Eq.  (6)  can  be  rephrased  as 

<t  =  M(D) :  &  (7) 

Clearly,  M(D)  is  the  transformation  tensor  relating  the  homogenized  and  effective  stresses.  More¬ 
over,  from  Eq.  (6),  we  obtain  the  elastic -damage  stiffness  tensor  C: 

C  =  M(D)  •  C°  (8) 

A  rational  choice  of  M(D)  may  be  M(D)  =  I  -  D ,  with  I  denoting  the  rank-four  unit  tensor.  If  this 
choice  is  made,  then  Eq.  (7)  and  (8)  can  be  recast  as  (see  also  Krajcinovic  (1989)) 

&  =  [I  -  D]-1  :  a  (9) 

C  =  [I  -  D]  •  C°  or  D  =  I-CC°-‘  (10) 

It  is  worth  noting  that  Eq.  (10)  shows  the  direct  relationship  between  the  damage  tensor  D  and  the 
elastic-damage  stiffness  tensor  C  .  Therefore,  damage  tensor  D  must  be  isotropic  (not  necessarily 
reducing  to  a  scalar  representation)  for  an  isotropic  damage  state.  Moreover,  the  appropriate  ex¬ 
pression  for  the  effective  stress  concept  is  Eq.  (7)  or  (9),  not  Eq.  (4)  even  in  the  case  of  isotropic 
damage.  On  the  other  hand,  D  and  hence  C  become  anisotropic  in  the  case  of  an  anisotropic  damage 


state. 
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In  this  section,  we  apply  some  results  from  micromechanical  analyses  for  two-dimensional 
microcrack-weakened  brittle  solids  under  tension  (Horii  and  Nemat-Nasser  (1983),  Sumarac  and 
Krajcinovic  (1987)).  Applications  can  be  made  to  brittle  matrix  fiber  or  particle  composites  as 
well  as  ceramics,  etc.  It  is  assumed  that  the  virgin  material  is  isotropic  and  linearly  elastic,  and 
that  line  microcracks  are  always  open  and  perfectly  randomly  distributed  in  all  directions  with 
uniform  size  2a.  The  microcrack  size  2a,  however,  may  grow  uniformly  during  the  loading  history. 
Therefore,  the  overall  material  response  remains  isotropic.  Furthermore,  the  overall  elastic-damage 
compliance  tensor  S  takes  the  form: 

S  =  S°  +  S*  (11) 

where  S°  is  the  virgin  elastic  compliance  (i.e.,  inverse  of  C°)  and  S*  denotes  the  damage-induced 
additional  inelastic  compliance. 

For  two-dimensional  problems,  it  is  more  convenient  to  use  the  Voigt’s  notation  and  express  the 
compliance  tensor  as  a  three  by  three  matrix.  Based  on  microcrack  opening  displacement  formulas, 
S*  can  be  micromechanically  computed  by  (assuming  plane  stress  condition): 

IN  a2  r 

s'  =  ~etJ0  mie  <12> 

where  N  =  the  number  of  microcracks  per  unit  cell  surface,  A  =  the  surface  area  of  representative 
unit  cell,  F  =  the  damaged  overall  “Young’s  modulus”,  and  g{9)  is  a  tensorial  function  of  the 
orientation  angle  6  ;  see  Sumarac  and  Krajcinovic  (1989,  Eq.  (17)).  Notice  that  S*  is  indeed  an 
isotropic  tensor. 


In  the  case  of  the  Taylor’s  model,  microcrack  interaction  is  completely  ignored.  By  actually 
evaluating  the  integral  in  Eq.  (12),  we  obtain  the  overall  elastic-damage  moduli  for  plane  stress  as 
follows  (see  Sumarac  and  Krajcinovic  (1989,  Eq.  (22))) 


E  v  _  1 

E0  va  1  +  u> 


(13) 


where  uj  =  Nna2/A  is  the  non-dimensional  microcrack  density  parameter,  and  E0  and  vD  are  the 
virgin  elastic  moduli.  Clearly,  u  <  v0  and  therefore  the  scalar  damage  representation  discussed 
in  Sec.  2  is  inadequate.  It  is  also  noteworthy  from  Eq.  (12)  that  S*  is  factually  a  function  of 
the  microcrack  density  parameter  u  .  Nevertheless,  one  should  not  confuse  the  scalar  microcrack 
density  parameter  u>  with  the  isotropic  damage  tensor  D  . 


Thus,  it  is  actually  possible  derive  the  fourth-order  damage  tensor  D  (or  D)  from  micromechan¬ 
ics  (i.e.,  the  microcrack-induced  inelastic  compliance  S').  Obviously,  if  microcracks  are  spatially 
randomly  distributed  in  a  representative  volume  element,  then  S’  and  hence  D  (or  D)  are  isotropic 
(given  linearly  elastic  and  isotropic  C°).  If  there  are  preferred  orientations  in  microcrack  distribu¬ 
tion,  then  S*  will  be  anisotropic  and  therefore  damage  tensor  D  (or  D)  becomes  anisotropic.  Finally, 
one  notes  that  if  S’  =  d  S°  then  D  =  d  I ,  and  hence  the  scalar  damage  model  discussed  in  Sec.  2  is 
recovered. 


t 
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1.5.  Conclusion 

I  The  present  work  analyzes  the  implication  and  limitation  of  some  scalar  damage  models.  It 

is  demonstrated  that  isotropic  damage  does  not  necessarily  imply  scalar  damage  representation  in 
general.  The  notion  of  isotropic  and  anisotropic  damage  variables  in  continuum  damage  mechan¬ 
ics  is  discussed.  In  addition,  some  results  from  micromechanical  analyses  are  applied  to  show  the 
*  direct  relationship  between  the  fourth-order  damage  tensor  D  (or  6)  and  the  damage-induced  com¬ 

pliance  tensor  S* .  It  is  concluded  that,  even  for  isotropic  damage,  one  should  employ  an  isotropic 
fourth-order  damage  tensor  (not  a  scalar  damage  variable)  to  characterize  the  state  of  damage  in 
I  materials.  In  general,  however,  a  damage  tensor  is  anisotropic.  Finally,  it  is  preferable  to  derive  gen¬ 

eral  anisotropic  damage  tensors  from  micromechanical  formulations,  instead  of  postulating  them 
arbitrarily. 

I 


» 


» 


» 


I 


I 


I 


On  the  Notion  of  Isotropic  and  Anisotropic  Damage  Variables 


9 


1.6.  References 

1.  Budiansky,  B.  and  R.  J.  O’Conell,  (1976),  “Elastic  Moduli  of  a  Cracked  Solid”,  Int. 

J.  Solids  &  Struct .,  Vol.  12,  pp.  81-97. 

2.  Chaboche,  J.  L.,  (1979),  “Le  Concept  de  Contrainte  Effective  Applique  a  L’elasticite  et 

a  la  Viscoplasticite  en  Presence  d’un  Endommagement  Anisotrope”,  Mechanical  Be¬ 
havior  of  Anisotropic  Solids,  Proc.  Euromech  Colloque  115,  June  1979,  Edited  by 
J.  P.  Boehler,  pp.  737-760,  Martinus  Nijhoff. 

3.  Cordebois,  J.  P.  and  F.  Sidoroff,  (1979),  “Damage  Induced  Elastic  Anisotropy”, 

Mechanical  Behavior  of  Anisotropic  Solids,  Proc.  Euromech  Colloque  115,  June  1979, 
Edited  by  J.  P.  Boehler,  pp.  761-774,  Martinus  Nijhoff. 

4.  Cordebois,  J.  P.  and  F.  Sidoroff,  (1982), “Endommagement  Anisotrope enElasticite 

et  Plasticite”,  J.  Mech.  Theor.  Appl.,  No.  special,  pp.  45-59. 

5.  Fanella,  D.  and  D.  Krajcinovic,  (1988),  “A  Micromechanical  Model  for  Concrete 

in  Compression”,  Eng.  Fract.  Mech.,  Vol.  29,  No.  1,  pp.  49-66. 

6.  Horii,  H.  and  S.  Nemat-Nasser,  (1983),  “Overall  Moduli  of  Solids  with  Microcracks: 

Load  Induced  Anisotropy”,  J.  Mech.  Phys.  Solids,  Vol.  31,  No.  2,  pp.  155-171. 

7.  Horii,  H.  and  S.  Nemat-Nasser,  (1985),  "Elastic  Fields  of  Interacting  Inhomogenei¬ 

ties”,  Int.  J.  Solids  &  Struct.,  Vol.  21,  No.  7,  pp.  731-745. 

8.  Ju,  J.  W.,  (1989a),  “On  Energy-Based  Coupled  Elastoplastic  Damage  Theories:  Consti¬ 

tutive  Modeling  and  Computational  Aspects”,  Int.  J.  Solids  &  Struct.,  Vol.  25,  No.  7, 
pp.  803-833. 

9.  Ju ,  J .  W.,  (1989b),  “On  Two-Dimensional  Self-Consistent  Micromechanical  Damage  Mod¬ 

els  for  Brittle  Solids”,  Int.  J.  Solids  &  Struct.,  accepted  for  publication,  in  press. 

10.  Kachanov,  L.  M.,  (1958),  “Time  of  the  Rupture  Process  under  Creep  Conditions”,  Izvestia 

Akademii  Nauk,  USSR,  Nauk  8,  pp.  26-31. 

11.  Kachanov,  L.  M.,  (1986),  “Introduction  to  Continuum  Damage  Mechanics”,  Martinus 

Nijhoff  Publishers,  Dordrecht,  Netherlands. 


> 


I 


I 


i 


» 


I 


t 


I 


p 


p 


On  the  Notion  of  Isotropic  and  Anisotropic  Damage  Variables  10 

12.  Kachanov,  M.,  (1980),  “A  Continuum  Model  of  Medium  with  Cracks”,  J.  Eng.  Mech. 

Div.,  ASCE,  Vol.  106,  pp.  1039-1051. 

13.  Kachanov,  M.,  (1987),  “Elastic  Solids  with  Many  Cracks:  A  Simple  Method  of  Analysis”, 

/nt.  /.  Solids  &  Struct.,  Vol.  23,  No.  1,  pp.  23-43. 

14.  Krajcinovic,  D.,  (1984),  “Continuum  Damage  Mechanics”,  Appl.  Mech.  Rev.,  Vol.  37, 

Jan.,  pp.  1-6. 

15.  Krajcinovic,  D.,  (1985),  “Constitutive  Theories  for  Solids  with  Defective  Microstruc¬ 

ture”,  in  Damage  Mechanics  and  Continuum  Modeling ,  ed.  N.  Stubbs  and  D.  Kra- 
jcinivic,  ASCE,  pp.  39-56. 

16.  Krajcinovic,  D.,  (1986),  “Update  to  Continuum  Damage  Mechanics”,  Appl.  Mech.  Up¬ 

date,  pp.  403-406. 

17.  Krajcinovic,  D.,  (1989),  “Damage  Mechanics”,  Mech.  of  Materials,  to  appear. 

18.  Krajcinovic,  D.  and  D.  Fanella,  (1986),  “A  Micromechanical  Damage  Model  for 

Concrete”,  Eng.  Fract.  Mech.,  Vol.  25,  No.  5/6,  pp.  585-596. 

19.  Krajcinovic,  D.  and  G.  U.  Fonseka,  (1981),  “The  Continuous  Damage  Theory  of 

Brittle  Materials”,  J.  Appl.  Mech.,  Vol.  48,  pp.  809-824. 

20.  Krajcinovic,  D.  and  D.  Sumarac,  (1989),  “A  Mesomechanical  Model  for  Brittle 

Deformation  Processes:  Part  I”,  J.  Appl.  Mech.,  Vol.  56,  No.  3,  pp.  51-56. 

21.  Lemaitre,  J.,  (1984),  “How  to  Use  Damage  Mechanics”,  Nuclear  Eng.  and  Design, 

Vol.  80,  pp.  233-245. 

22.  Lemaitre,  J.  and  J.  L.  Chaboche,  (1978),  “Aspect  Phenomenologiquede  la  Rupture 

par  Endommagement”,  J.  de  Mech.  Applique,  Vol.  2,  pp.  317-365. 

23.  Lemaitre,  J.  and  J.  L.  Chaboche,  (1985), Mechanique des Materiaux Solids,Dunod, 

Paris. 

24.  Mazars,  J.,  (1984),  “Application  de  la  Mechanique  de  1 ’Endommagement  au  Comporte- 

ment  non  Lineaire  et  a  la  Rupture  du  Beton  de  Structure”,  These  de  Doctorate  d’Etat, 
L.M.T.,  Universite  Paris,  France. 


> 


On  the  Notion  of  Isotropic  and  Anisotropic  Damage  Variables 


11 


25.  Murakami,  S.  and  N.  Ohno,  (1981),  “A  Continuum  Theory  of  Creep  and  Creep  Dam¬ 

age”,  in  Creep  of  Structures,  IUTAM  Symp.,  Ed.  A.R.S.  Ponter,  Springer- Verlag,  Berlin, 
pp.  422-444. 

26.  Rabotnov,  I.  N.,  (1963),  “On  the  Equations  of  State  for  Creep”,  in  Progress  in  Applied 

Mechanics-The  Prager  Anniversary  Volume ,  MacMillan,  New  York,  pp.  307-315. 

27.  Sumarac,  D.  and  D.  Krajcinovic,  (1987),  “A  Self-Consistent  Model  for  Microcrack- 

Weakened  Solids”,  Mech.  Mater.,  Vol.  6,  pp.  39-52. 

28.  Sumarac,  D.  anD.  Krajcinovic,  (1989),  “A  Mesomechanical  Model  for  Brittle  De¬ 

formation  Processes:  Part  IT’,  J.  Appl.  Mech.,  Vol.  56,  No.  3,  pp.  57-62. 

29.  Vakulenko,  A.  A.  and  M.  Kachanov,  (1971),  “Continuum  Theory  of  Cracked  Me¬ 

dia”,  Izv.  ANSSSR,  Mekh.  Tverdogo  Tela,  Vol.  4,  pp.  159-166. 


12- 

PART  II: 

A  Micromechanical  Damage  Model  for  Uniaxially  Reinforced 
Composites  with  Interfacial  Arc  Microcracks 


H.O.  Abstract 

A  micromechanical  damage  model  is  presented  in  this  work  for  fiber/matrix  interfacial  micro¬ 
crack-weakened  uniaxially  reinforced  brittle  matrix  composites.  It  is  assumed  that  the  undamaged 
fibers  and  matrix  materials  are  linearly  elastic  and  isotropic.  All  microcracks  are  assumed  to  occur 
in  the  fiber/matrix  interfaces,  and  are  modeled  as  arc  microcracks  under  “cleavage  1”  deforma¬ 
tion  processes.  Microcrack-induced  inelastic  strains  and  overall  elastic-damage  compliances  are 
analytically  derived  based  on  micromechanical  bi-material  (interfacial)  arc-microcrack  opening 
displacements  and  probabilistic  distributions  of  arc  microcracks.  In  addition,  thermodynamic  ba¬ 
sis  is  rendered.  The  present  work  is  entirely  different  from  existing  phenomenological  continuum 
damage  models  for  composite  materials  since  it  is  based  on  micromechanical  derivations.  Both 
“ stationary ”  and  “ process "  damage  models  are  given  in  the  present  work.  In  particular,  microcrack 
kinetic  equations  are  constructed  based  on  micromechanical  fracture  criterion  and  microstructural 
geometry  in  a  representative  volume  element.  Simple  and  efficient  computational  algorithms  as 
well  as  some  numerical  uniaxial  tension  tests  are  also  presented  to  illustrate  the  proposed  microme¬ 
chanical  damage  model  for  composites.  It  is  emphasized  that  the  present  work  does  not  employ 
any  arbitrary  (fitted)  “material  constants”. 
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II.l.  Introduction 

Initiated  by  Kachanov  (1958)  and  Rabotnov  (1963)  for  one-dimensional  creep  damage  of  met¬ 
als,  continuum  damage  mechanics  has  been  extensively  explored  and  applied  by  many  outstanding 
applied  mechanics  researchers.  Basically,  there  are  two  types  of  continuum  damage  formulations  - 
phenomenological  and  micromechanical  damage  models.  Most  existing  work  are  classified  as  phe¬ 
nomenological  damage  models;  see,  e.g.,  Krajcinovic  (1984,1986,1989)  and  Bazant  (1986)  for  a 
comprehensive  literature  review.  There  are,  however,  some  micromechanical  “stationary”  or  “pro¬ 
cess”  damage  models  proposed  in  the  literature;  see,  e.g.,  Budiansky  and  O’Connel  (1976),  Horii 
and  Nemat-Nasser  (1983,  1985),  Kachanov  (1987),  Krajcinovic  and  Fanella  (1986),  Fanella  and 
Krajcinovic  (1988),  Sumarac  and  Krajcinovic  (1987),  and  Ju  (1989b). 

In  particular,  excellent  studies  on  damage  mechanics  in  modem  composite  materials  were  pre¬ 
sented  by  Weitsman  (1987a,1987b,1988),  Talreja  (1985a, 1985b, 1985c, 1986,1987),  Allen  et  al. 
(1987a,1987b,1987c),  and  Harris  et  al.  (1987)  for  distributed  matrix  cracks  and  delamination 
within  the  framework  of  phenomenological  damage  models.  On  the  other  hand,  some  noteworthy 
micromechanical  (primarily  “stationary”)  damage  models  for  composites  were  proposed  by,  e.g., 
Wang  et  al.  (1984),  Laws  et  al.  (1983)  and  Hashin  (1985)  for  transverse  (parallel)  matrix  crack 
systems.  It  is  noted  that  existing  phenomenological  continuum  damage  models  of  Weitsman,  Tal¬ 
reja  and  Allen  et  al.  employed  either  vector-valued  or  second-rank  (symmetric  or  non-symmetric) 
“damage  tensors”  (treated  as  internal  state  variables)  to  characterize  the  state  of  damage  in  compos¬ 
ite  materials.  However,  a  vector  or  a  second-rank  damage  tensor  is  inherently  incapable  of  describ¬ 
ing  general  anisotropy  in  composites.  An  appropriate  description  of  anisotropic  damage  generally 
involves  &  fourth-rank  (or  even  eighth-rank)  damage  tensor  representation;  see,  Chaboche  (1979), 
Cordebois  and  Sidoroff  (1979,1982),  Ju  (1989a),  and  Krajcinovic  (1989)  for  further  remarks.  In 
addition,  in  spite  of  the  solid  and  attractive  thermodynamic  basis,  the  specific  functional  forms  of 
the  Helmholtz  or  Gibbs  free  energy  potentials  in  phenomenological  damage  models  are  to  some 
extent  arbitrary  (heuristic).  Therefore,  the  resulting  overall  stiffness-damage  relationships  and 
stress-strain  laws  are  also  somewhat  arbitrary.  Moreover,  in  order  to  have  constitutive  predictive 
capability,  phenomenological  damage  models  must  arbitrarily  or  empirically  postulate  functional 
forms  for  damage  “evolution  (growth)  equations”.  In  the  thermodynamic  free  energy  potentials 
and  damage  evolution  equations,  nevertheless,  existing  phenomenological  damage  models  typi¬ 
cally  rely  on  the  use  of  many  (perhaps  up  to  100  or  even  200)  fitted  “material  constants”.  As  a 
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consequence,  it  becomes  very  difficult  to  identify  these  fitted  constants  from  actual  experimental 
data  of  composites. 

Hence,  as  pointed  out  by  Krajcinovic  and  Fanella  (1986)  and  Weitsman  (1988),  micromechan¬ 
ical  damage  theories  for  composites  are  warranted  to  incorporate  microstructural  geometry,  mi¬ 
cromechanical  deformations  and  microcrack  growth  into  the  damage  mechanics  framework.  Most 
currently  available  micromechanical  damage  models  for  composites  focused  on  the  effects  of  trans¬ 
verse  stationary  matrix  cracks  on  overall  (damaged)  compliance  tensors.  The  present  paper,  on  the 
other  hand,  considers  the  damage  effects  on  uniaxially  reinforced  composites  due  to  the  existence 
and  growth  of  an  ensemble  of  randomly  oriented  microcracks  at  the  fiber/matrix  interfaces  within 
the  context  of  micromechanical  damage  mechanics.  This  circumstance  corresponds  to  an  ensem¬ 
ble  of  randomly  distributed  arc  microcracks  at  the  interfaces  between  cyclindrical  inclusions  and 
extended  exterior  regions  under  plane  strain.  In  the  case  of  a  single  arc  crack  at  the  fiber/matrix 
interface  subjected  to  remote  tension  field  (plane  strain),  solutions  are  available  in  England  (1966), 
Perlman  and  Sih  (1967),  Toya  (1974),  and  Piva  (1982).  In  particular,  Toya’s  solution  is  the  most 
suitable  for  our  damage  mechanics  formulation  because  it  provides  analytical  (micromechanical) 
expressions  for  microcrack  opening  displacements  and  arc  microcrack  fracture  energy  criteria;  see 
also  remarks  made  in  Weitsman  (1988). 

In  the  present  work,  it  is  assumed  that  the  undamaged  fibers  and  matrix  materials  have  different 
properties  but  are  both  linearly  elastic  and  isotropic.  All  arc  microcracks  are  assumed  to  occur  in 
the  fiber/matrix  interfaces  under  “cleavage  1”  deformation  processes  (i.e.,  due  to  pre-existing  arc 
microcracks).  New  microcrack  nucleations  (“cleavage  2”  processes)  are  not  considered  in  this 
work.  Moreover,  at  this  stage  of  the  development,  arc-microcrack  interactions  are  assumed  to  be 
negligible  (Taylor’s  model);  i.e.,  dilute  microcrack  concentrations  are  considered.  In  the  absence  of 
a  micromechanical  solution  for  partially  closed  or  entirely  closed  arc  microcrack  (namely,  a  “mixed 
crack-and-contact  problem”),  we  assume  that  all  arc  microcracks  are  entirely  open  under  a  remote 
tension  field. 

An  outline  of  this  paper  is  as  follows.  The  representation  of  the  fourth-rank  damage  and  ther¬ 
modynamic  basis  for  microcrack-weakened  brittle  composites  are  presented  in  Sec.  n.2  within 
the  context  of  homogenization  for  multiphase  inhomogeneous  materials.  It  is  assumed  that  dis¬ 
tributed  arc  microcrack  concentration  justifies  the  use  of  effective  continuum  medium  theory.  Based 
on  Toya’s  (1974)  micromechanical  solution  of  bi-material  (interface)  arc-crack  opening  displace¬ 
ments,  damage-induced  inelastic  strains  and  inelastic  compliances  are  systematically  derived  in 
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Sec.  11.3  for  an  ensemble  of  randomly  oriented  open  arc  microcracks.  This  corresponds  to  a  sta¬ 
tionary  (non-process)  micromechanical  damage  model.  In  Sec.  H.4,  microcrack  growth  (evolution) 
is  considered  based  on  Toya’s  (1974)  micromechanical  fracture  criterion  for  a  single  arc  microc¬ 
rack  under  uniaxial  tension.  The  extension  to  account  for  biaxial  tension  loadings  can  be  readily 
made.  “Stable”  and  “unstable”  domains  are  identified  for  stationary  and  propagating  arc  micro¬ 
cracks,  respectively.  As  a  consequence,  a  “process  model”  is  rendered.  It  is  emphasized  that  the 
present  work  does  not  employ  any  arbitrary  (fitted)  “material  constants”.  Simple  and  efficient  com¬ 
putational  algorithms  are  given  in  Sec.  n.5.  In  addition,  some  numerical  uniaxial  tension  tests  are 
presented  in  Sec.  H.5  to  illustrate  the  potential  capability  of  the  proposed  micromechanical  damage 
model  for  composites. 
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II.2.  Thermodynamic  basis 

In  the  case  of  general  anisotropy,  a  scalar,  a  vectorial,  or  a  second-rank  tensorial  representation 
of  damage  variable  is  not  adequate;  see  Krajcinovic  (1989)  for  more  details.  In  this  work,  we 
employ  a  fourth-rank  anisotropic  damage  tensor  to  represent  the  state  of  damage  in  anisotropic 
composite  materials.  It  is  worth  mentioning  that  the  fourth-rank  damage  tensor  employed  has  an 
appealing  correspondence  with  the  fourth-rank  overall  stiffness  (or  compliance)  tensor. 

Within  the  framework  of  homogenization  concept  for  inhomogeneous  effective  continuum 
medium,  one  may  define  the  homogenized  Gibbs  free  energy  as 

X  =  i<r:[S0(I  +  D)]:<7  (1) 

where  er  is  the  volume-average  stress  tensor  (Hill  (1965)),  S°  is  the  undamaged  linear  elastic  com¬ 
pliance  (obtained  by  the  rule  of  mixture)  of  a  composite  material,  I  is  the  fourth-rank  unit  tensor,  D 
denotes  the  fourth-rank  damage  tensor,  and  “  :  ”  denotes  the  tensor  contraction  operation.  By  the 
Clausius-Duhem  inequality  for  isothermal  process,  we  have  (with  e  denoting  the  volume-average 
strain) 

X  -  dr :  e  >  0  (2) 

The  standard  Coleman’s  method  then  leads  to  the  following  macroscopic  stress-strain  law  and  over¬ 
all  elastic-damage  compliance  tensor  S  : 

c  =  [S°  •  (I  +  D)] :  cr  (3) 

S  =  S°-(I  +  D)  (4) 

together  with  the  damage  dissipative  inequality: 

(5) 

From  Eq.  (5),  it  is  observed  that  the  evolution  Js  plays  an  essential  role  in  microcrack  energy 
dissipation  and  growth.  Further,  Eq.  (3)  is  entirely  consistent  with  the  “effective  stress”  concept 
in  continuum  damage  mechanics  (Kachanov  (1958),  Rabotnov  (1963),  Lemaitre  and  Chaboche 
(1978)).  During  a  damage  loading  process,  the  total  strain  tensor  e  is  amenable  to  an  additive 
decomposition: 


€  =  e'  +  e 


(6) 
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where  ee  and  e *  denote  the  elastic  and  inelastic  (damage-induced)  strains,  respectively.  In  addition, 
it  is  assumed  that  em  «  0  upon  complete  unloading;  namely,  the  residual  strain  at  zero  stress  is 
negligible  for  brittle  composite  materials.  The  elastic-damage  compliance  tensor  is  also  suitable 
for  an  additive  decomposition: 

S  =  S°  +  S*  (7) 

where  S*  signifies  the  damage-induced  inelastic  compliance  (Mura  (1982),  Horii  and  Nemat-Nasser 
(1983)).  In  fact,  from  Eqs.  (4)  and  (7),  the  relationship  between  S*  and  D  can  be  formally  expressed 
as 

S*  =  S°  •  D  (8) 

Therefore,  if  one  can  micromechanically  derive  the  damage-induced  inelastic  compliance  S*, 
then  one  can  explicitly  express  the  fourth-rank  damage  tensor  D  by  means  of  micromechanics 
in  stead  of  postulating  D  and  6  heuristically.  This  is  precisely  the  motivation  of  the  following 
development. 
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IU.  Microcrack-induced  inelastic  strains  and  overall  compliances 

In  this  section,  damage-induced  inelastic  strains  and  inelastic  compliances  are  derived  for  an 
ensemble  of  randomly  distributed  iiber/matrix  interfacial  arc  microcracks.  Microcrack  interactions 
are  neglected  at  this  stage  of  our  development,  and  shall  be  subject  of  future  study.  According,  only 
“Taylor’s  model”  (not  self-consistent  model)  is  constructed.  The  microcrack  opening  displacement 
formulas  for  a  single  interfacial  arc  microcrack  are  based  on  Toya’s  (1974)  solution.  All  interfacial 
microcracks  are  assumed  to  be  entirely  open.  If  micromechanical  solutions  for  partially  closed 
or  entirely  closed  interfacial  arc  microcrack  become  available  in  the  future,  however,  the  present 
framework  can  be  readily  modified  to  accommodate  those  circumstances.  The  development  in  this 
section  can  be  classified  as  “stationary”  (or  “non-process")  model  according  to  the  terminology  of 
Krajcinovic  and  Fanella  (1986). 

II.3.1.  Microcrack-induced  inelastic  strains 

In  Toya’s  (1974)  solution,  the  fiber  (inclusion)  and  the  matrix  materials  are  assumed  to  have 
different  linear  elastic  properties.  Nevertheless,  both  fiber  and  matrix  constituents  are  assumed  to 
be  homogeneous  and  isotropic.  Toya  (1974)  actually  provides  solutions  for  stresses,  displacements, 
and  debonding  criteria  for  an  open  arc  crack  at  the  bi-material  interface  under  remote  uniaxial  and 
biaxial  tension  loadings.  In  what  follows,  for  simplicity,  we  only  consider  the  case  of  uniaxial 
tension  loadings. 

As  remarked  by  Toya  (1974),  both  stresses  and  displacements  oscillate  violently  at  the  imme¬ 
diate  regions  near  the  crack  tips.  This  is  quite  typical  for  the  mixed  boundary-value  problem  for 
interfacial  cracks  at  bi-material  boundaries.  However,  the  extent  of  the  oscillating  regions  is  very 
small.  Consequently,  Toya  (1974)  concludes  that  his  solutions  provide  a  good  approximation  to 
the  physical  state  of  the  body  at  the  bi-material  interface  except  in  the  immediate  vicinity  of  the 
crack  tips.  Although  some  recent  results  of  oscillation-free  bi-material  stress  and  displacement 
analyses  were  reported  in  several  technical  symposiums  and  journals  (see,  e.g.,  Qu  and  Bassani 
(1989a,  1989b),  Bassani  and  Qu  (1989)),  they  were  derived  for  line  (straight)  microcracks  at  the 
interface  of  two  dissimilar  semi-infinite  materials,  not  valid  in  our  arc  crack  problem.  Therefore, 
Toya’s  solution  is  adopted  in  the  present  work. 

The  global  (unprimed)  and  local  (primed)  Cartesian  coordinate  systems  as  well  as  the  local 
polar  coordinate  system  (at  a  typical  arc  point)  are  shown  in  Figure  1.  In  particular,  a  denotes  the 
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half-angle  expanded  by  an  arc  microcrack,  a  denotes  the  radius  of  the  fiber  (inclusion),  <j>  denotes 
the  angle  between  the  x'-axis  and  y-axis,  and  ip  (=  tt/2  -  <p)  signifies  the  angle  between  the  y-axis 
and  y'-axis.  The  uniaxial  tension  p  is  applied  in  the  y-axis  direction.  Counterclockwise  direction 
is  taken  as  positive,  and  0  is  measured  from  the  x'-axis.  In  addition,  n'  and  (u'r ,  u'e)  represent  the 
outward  unit  normal  vector  and  the  polar  coordinates  at  a  typical  point  along  the  arc,  respectively. 
The  expressions  for  u'T  and  u'e  under  remote  uniaxial  tension  are  given  by  Eq.  (3.57)  in  Toya 
(1974)  by  means  of  a  complex  variable  form.  Note  that  the  real  part  is  u'T  and  the  imaginary  part 
is  u'g .  Let  us  first  transform  the  local  polar  coordinates  at  a  typical  arc  point  to  the  local  Cartesian 
coordinate  defined  at  the  mid-point  of  the  arc: 


=  u'r  cos  9  —  u'e  sin  9 

(9a) 

=  u'r  sin  9  +  u'e  cos  6 

(9b) 

Furthermore,  the  Cartesian  components  of  the  unit  outward  normal  vector  at  a  typical  arc  point  are: 

n'x=cos0  ;  n'x  =  cos  0  (10) 

Since  the  present  work  is  concerned  with  the  mechanical  behavior  of  interfacial  microcrack- 
weakened  uniaxially  reinforced  brittle  composites,  it  is  appropriate  to  focus  on  the  c  ise  of  plane 
strain  (orthogonal  to  the  fiber  direction).  The  virgin  brittle  composite  material  under  consider¬ 
ation  is  actually  transversely  isotropic.  Thus,  in  a  plane  strain  framework,  the  virgin  composite 
material  is  isotropic.  However,  if  the  composite  contains  an  ensemble  of  microcracks,  it  may  be¬ 
come  anisotropic,  depending  on  the  microcrack  sizes  and  orientations.  By  the  standard  definition 
of  Cartesian  strain  tensor,  the  inelastic  plane  strain  components  (in  Voigt’s  notation)  due  to  a  single 
(fc-th)  arc  microcrack  take  the  form: 


t?r  =  c-j>r  =  ljyxmsl,d0'.?.j 

r° 

’  (u'T  cos2  6  -  u'e  sin  0  cos  0)  d9 

— a 

(11) 

efy  =  c«>y  =  jJ"U',sineiO  =  jj 

f  (u'r  sin2  9  +  u'e  sin  0  cos  9)  d9 

—a 

(12) 

c*(*)'  s  2e*(*),  f  (u'x  sin  6  +  u'y  cos  9)d9  - 

^  J —a  ^ < 

/  (u'r  sin  20  —  u'e  cos  20)  d0 

J  —  Or 

(13) 

where  A  is  the  surface  area  of  a  representative  volume  element  in  two-dimension. 

In  the  above  equations,  it  is  implicitly  assumed  that  the  fc-th  arc  microcrack  is  entirely  open. 
That  is,  the  radial  crack  opening  displacement  u'r  is  always  positive.  Therefore,  there  exist  some 
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restrictions  on  the  arc  microcrack  size  2a  and  the  orientation  <f>;  see  Toya  (1974)  for  further  remarks. 
For  example,  in  the  case  of  the  epoxy  matrix  (shear  modulus  p\  =  346  KSI  or  2.39  GN/m2,  Poisson’s 
ratio  rji  =  0.35)  and  glass  fiber  (shear  modulus  p2  =6410  KSI  or  44.2  GN/m2,  Poisson’s  ratio 
Vi  =  0.22)  composite  material,  the  range  of  “entirely  open”  arc  microcracks  is  approximately 
defined  by: 

|0|+a<65°  (14) 

See  Fig.  2  in  Toya  (1974)  for  more  information  regarding  allowable  (<p,  a)  region. 

The  inelastic  strains  due  to  an  ensemble  of  non-interacting  arc  microcracks  can  be  evaluated 
by  performing  the  following  integration: 

e?  =  N  [  gflf1  W,a)  dl>  (15) 

Jn 

where  N  is  the  number  of  open  (active)  arc  microcracks;  i,  j  =  1, 2, 6;  P(0,  a)  is  a  joint  probability 
density  function  of  arc  microcrack  orientations  and  sizes;  Q  is  the  domain  of  all  open  (active) 
microcracks,  and  gW  is  the  component  of  the  following  local-global  transformation  matrix  (Horii 
and  Nemat-Nasser  (1983)) 

cos2  ip  sin2  ip  sin  2  ip 

[g(fc)]  =  sin 2  ip  cos2  ip  -sin  2V>  (16) 

-  5  sin  ftp  |  sin  ftp  cos  ftp 

Eq.  (15)  can  be  numerically  computed  by  double  Gauss  quadrature,  as  will  be  discussed  in  Sec. 
II.5. 

The  total  strain  components  can  be  obtained  by  adding  e*  to  the  elastic  contributions  t\ ,  with 
e*  expressed  as 

A  =  07) 

Remark  3.1.  In  the  case  of  biaxial  tension  loadings,  e*  can  be  obtained  by  the  same  procedure 
as  above.  The  only  modification  required  is  to  use  Eq.  (3.43)  (not  Eq.  (3.57))  in  Toya  (1974)  when 
integrating  arc  crack  opening  displacements.  | 


II 3.2.  Overall  elastic-damage  compliances 

%  To  derive  the  damage-induced  inelastic  compliance  matrix  S* ,  Eqs.  (1 1)— (13)  must  be  mod¬ 

ified.  The  key  step  is  to  construct  the  inelastic  strain-stress  relationship  (e*w'  vs.  a)  in  local 


t 
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Cartesian  coordinates.  However,  this  local  relationship  is  not  readily  available  from  (3.57)  in  Toya 
(1974)  since  crack  opening  displacements  are  related  to  global  tension  stress  p  instead  of  local 
stress  field  in  Toya’s  work.  The  local  stresses  can  be  easily  obtained  as  follows 


o\  =  oxx  =  p  COS2  4> 

(18) 

o'l  =  a'yy=P  sin2  <t> 

(19) 

'  '  p  • 

=  °xy  =  2  Sm2(^ 

(20) 

but  straightforward  derivation,  Eqs.  (1 1>— (13)  (due  to  a  single  arc  microcrack 

as 

|  ^2  \^J  (u'r« c°s2  ^  —  u'e,  si°  9  cos  0 )  d0  o\  | 

(21) 

C^k)  =  ^  [/  ( u'n  sin2  0  +  u'e,  sin  0  cos  0)  dO  o'  | 

(22) 

c^(fc) = ^  |  [y  2^  ~  u,«<  c°s 2^)  m  o\  | 

(23) 

where  the  summation  is  for  i  =  1, 2, 6,  and  the  expressions  for  u'ri  and  u'e,  are 

I 

uri  =  —  Ai<i  [sin  |(a  —  0)  sin  j(a  +  0)]  2  eA(,r_0,)  x 

^  1 1 1  —  (cos a  +  2A sin a)e2A(,r_0,)  +  (1  —  ifc)(l  +  4A2) sin2 a  1 


2  —  k  —  A:(cos  a  +  2 A  sin  a)e2A(,r-o) 


_  ~  fr)  2A(a-iQ 


cos 


2(1  -  k)^Ha_r) 


J  Uini(a+0)Jy 

V2  [smA(a  +  0)J/J  J 


sin  0  sin 


u'«i  =  —A\a  [sin  j(a  —  0)  sin  j(a  +  0)] 2  eA(*  o)  x 


—  (cos  a  +  2A  sin  a)e2A(,r-o)  +  (1  —  Jk)(l  +  4A2)  sin2  a  1 
2  —  k  —  k( cos  a  +  2A  sin  a)e2A(’r~°,)  k 


_  2(!  -  *0  2A(o-ir) 


2(1  -  k)  _2X(a-n) 


coS»]sin(l9-Aln[^7|])  + 

<’cos(i#-AinEfc|])]} 


(24) 


(25) 
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u'r 2  =  -A\a  [sin  j(q  -  0)  sin  5(0  +  0)] 5  ex^~Q)  x 

1  —  (cos  a  +  2A  sin  a)e2x^~a)  —  (1  —  Jfc)(l  +  4A2)  sin2  a 


-f[! 


2  —  k  —  k( cos  a  +  2A  sinaje2^"-0') 


1 

k  + 


2(1  -  k) 


e2A(a-ir)  sjn , 


f =0  -  A  In 

sin  5(0  -  0)’ 

sin  j(a  +  0) . 
‘sin  5(0  -  0)' 

V2 

f  50  —  A  In 

sin  5(0  +  0) . 

V 

u'n  =  -Aid  [sin  i(a  -  0 )  sin  A(a  +  0)] 5  eA(*_a)  x 

^  |  |1  —  (cos  a  +  2A  sin  ar)e2A(lr-o)  —  (1  —  fc)(l  4A2)  sin2  a 


2  —  k  —  k( cos  a  +  2A  sin  a)e2A^-“) 


1 

“  k  + 


+  2a_i)e„<„-„cos#]  sin  ^  _  A1„  - 

in*cos(^-A.ng|^|])]} 


2(1  -  k)  jmc-*) 


u'T 6  =  ~2A\a  [sin  A(o  —  0)  sin  5(0  +  0)P  eX{*  o)  x 


-I 


2(1  -  k) 


sin  0  cos 


\2  sin  ±(a  +  0)  J  / 


x  sin 


(l-t)(l+4*Wo  _  2(1  -  A)  - 

fc[  1  +  (cos  a  +  2A  sin  a)e2A(*  "°>]  k 

\  sin  i(a  +  0)  J  /  J 


u'b(>  =  -2Aia  [sin  |(a  —  0)  sin  j(a  +  0)] 2  eA(,r~o)  x 

f 2^  _ 

1  k  v2  lsin|(a  +  0)Jy 

_  [  (1  —  fc)(l  +  4A2)  sin2  a  _  2(1  -  *)»,„..)  ■ 

[fc[l +(cosa +  2Asina)e2A(,r-0>]  fc 

(f—IStel)} 


x  cos 


(26) 


(27) 


(28) 


(29) 


In  the  above  equations,  k ,  A,  and  Ai  are  given  material  properties  (related  to  elastic  shear  moduli 
and  Poisson’s  ratios  of  fiber  and  matrix),  and  are  defined  in  (3.9),  (3.3),  (3.39)  in  Toya  (1974).  For 
convenience,  they  are  given  in  what  follows: 


k  = 


_  .  ^  _  fll(l  +«2) 

1+1/  ’  _(/ll+«l/l2) 

Ki  =  3  -  4  T?1  ;  K2 


_  (P2  +  *2Pl) 
~  (Pi  +  M/^) 

3  —  4  772 


(30) 


(31) 


I 
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I 
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\ 
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i 


i 


A\  = 


k 

4 


A  = 


In  v 
2  7T 


1  +(C] 


01 


1  +  *2 

. 


Therefore,  Eqs.  (21H23)  can  be  rewritten  as  (t,  j  =  1, 2, 6) 


(32) 

(33) 


(34) 


where  S*fk^  is  the  local  two-dimensional  microcrack-induced  inelastic  compliance  component.  In 
particular,  the  inelastic  compliance  matrix  components  take  the  form: 


'  CL  fQ 

'-fc)  =  —  J  (u'ri  cos2  6  -  u'6i  sin  6  cos  9)  d9 

(35) 

i  CL 

S^-k)  =  —  j  (u'ri  sin2  9  +  u'ei  sin  9  cos  9)  d9 

(36) 

Sll k)  -  ^  J  (u'ri  sin  29  -  u'6i  cos  29)  d9 

(37) 

It  is  noticed  that  the  inelastic  compliance  matrix  S*^  due  to  a  single  arc  microcrack  is  in  fact 
non-symmetric\  i.e.,  S-jky  i  S‘,(ky  for  i  i  j.  The  integrals  in  Eqs.  (35)-(37)  can  be  effectively 
evaluated  by  means  of  numerical  integration  procedures;  see  Sec.  D.5  for  details.  In  addition,  in 
terms  of  global  Cartesian  coordinates,  S*}k)'  can  be  rephrased  as 


(38) 


Consequently,  the  total  inelastic  compliance  S*  due  to  an  ensemble  of  randomly  oriented  interfacial 
arc  microcracks  can  be  expressed  as 

=  N  f  S*»P(+,C)  to  (39) 

Jn 

From  Eqs.  (35)— (37),  it  is  observed  that  S*  is  actually  a  tensorial  function  of  the  mean  “arc  mi¬ 
crocrack  density  parameter”  <  u  >  ,  with  u  =  N(aa)2/A  .  Finally,  the  overall  elastic-damage 
compliance  S  is  obtained  simply  by  writing 


S  =  S°  +  S* 


(40) 


Therefore,  the  inelastic  compliances  and  fourth-rank  damage  tensors  of  interfacial  microcrack- 
weakened  brittle  composites  can  actually  be  derived  from  micromechanics,  instead  of  postulating 
them  arbitrarily  or  heuristically. 


i 
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It  is  characteristic  that  brittle  composite  solids  have  some  pre-existing  arc  microcracks  along 
the  fiber/matrix  interfaces  even  before  specimens  are  first  loaded.  Some  of  these  initial  microcracks 
may  grow  (propagate)  later  upon  application  of  adequate  loads.  In  this  section,  we  transform  the 
stationary  (“non-process”)  damage  model  presented  in  Sec.  II. 3  into  a  “process”  damage  model  so 
that  the  model  possesses  a  constitutive  predictive  capability.  That  is,  “cleavage  1”  (pre-existing) 
microcrack  growth  and  the  evolution  of  S  (or  S*)  are  considered  based  on  micromechanical  con¬ 
siderations.  The  importance  of  micromechanically-based  microcrack  growth  relations  were  also 
recognized  by  Weitsman  (1988).  It  is  recalled  that  although  some  existing  phenomenological  dam¬ 
age  models  are  in  principle  capable  of  providing  damage  evolution  equations  for  composites,  they 
may  require  some  100  (or  even  more)  fitted  “material  constants”  in  the  microcrack  evolution  laws. 

The  present  work  is  nor  restricted  to  monotonically  increasing  loads.  In  fact,  loading/unloading 
sequences  can  be  easily  accommodated  by  computing  and  checking  whether  there  is  undergoing 
microcrack  growth  (excluding  those  previously  propagating  and  currently  arrested  microcracks). 
If  there  is  no  (<f>,  a)  region  in  which  additional  microcrack  growth  is  now  taking  place,  then  the 
current  incremental  load  step  is  in  an  unloading  state.  Therefore,  “active  microcrack  growth”  is 
the  valid  current  loading  condition,  regardless  of  prior  existence  of  certain  {<j>,  a)  regions  where 
microcracks  previously  experienced  growth.  Moreover,  the  damage-induced  inelastic  compliance 
S*  takes  the  ftrm: 

s*-s;+s;  +  s}  (4i) 

where  S*  =  the  compliance  contribution  from  undergoing  microcrack  growth,  S*  =  the  contribution 
from  stationary  microcracks  having  initial  sizes,  and  SJ  =  the  contribution  from  arrested  microcrack 
growth  due  to  previous  microcrack  growth.  In  particular,  if  S*  =  0,  then  the  current  load  level  p 
is  not  high  enough  to  cause  further  damage,  and  therefore  all  existing  microcracks  are  arrested. 
Finally,  S*  is  added  to  S°  to  obtain  the  overall  compliance  S. 

n.4.1.  Interfacial  microcrack  fracture  criterion 

The  mixed  mode  bi-material  fracture  criterion  was  provided  in  Eq.  (4.8)  in  Toya  (1974)  for  a 
non-interacting,  entirely  open  arc  microcrack  along  the  fiber/matrix  interface  under  uniaxial  ten¬ 
sion.  It  is  implicitly  assumed  that  the  bonding  strength  between  the  fiber  and  matrix  is  sufficiently 
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small  compared  with  the  fracture  toughness  of  the  matrix,  so  that  an  initial  (existing)  arc  micro¬ 
crack  will  always  grow  along  the  bi-material  interface.  Thus,  microcrack  kinking  into  the  matrix  is 
not  considered  here.  Toya’s  bi-material,  local  stress  fracture  criterion  for  a  tip  of  an  arc  microcrack 
under  uniaxial  tension  reads: 

\p2kaAx(\  +  4\2)kN0N0  sin  a  e2X^~a)  =  27 iJ  (42) 

16 

where  7'^  is  the  specific  surface  tension  energy  of  the  interface  (i.e.,  critical  energy  release  rate),  and 
No  together  with  its  complex  conjugate  No  are  functions  of  4>,a,  and  elastic  material  properties  (see 
Eq.  (4.9)  in  Toya  (1974)).  Eq.  (42)  is  only  valid  for  an  open  microcrack;  i.e.,  u'r  >  0  is  required 
for  any  9  value  along  an  arc  microcrack.  In  addition,  one  tip  of  a  microcrack  may  reach  the  fracture 
criterion  (4.2)  before  the  other  tip  does.  Hence,  one  tip  may  be  propagating  while  the  other  tip  is 
arrested. 

When  the  energy  release  rate  reaches  its  critical  value  (or  when  tension  reaches  a  critical  value 
Per),  an  arc  microcrack  may  grow  in  a  stable  or  unstable  fashion,  according  to  Fig.  3  in  Toya 
(1974).  Within  a  limited  range  of  the  ( <f> ,  a)  region,  Eq.  (42)  may  have  two  solutions  for  a  and 
4>.  Thus,  one  tip  or  both  tips  may  actually  grow  in  a  stable  manner.  This  implies  that  the  final 
microcrack  size  qj  and  orientation  <f>{  can  be  analytically  obtained.  However,  for  other  ranges  of 
the  (4>,  a)  region,  there  is  only  cne  solution  for  Eq.  (42).  Therefore,  an  arc  microcrack  may  grow 
in  an  unstable  manner  outside  the  limited  range  of  the  two-solution  ( <f> ,  a)  domain.  Moreover,  with 
further  increase  of  the  tensile  stress  p  (>  p^),  even  an  originally  stable  arc  microcrack  will  grow 
continuously  in  an  unstable  manner. 

In  the  case  of  many  randomly  distributed  open  arc  microcracks,  Eq.  (42)  is  systematically 
checked  for  every  permissible  microcrack  orientation  <f>  and  size  a  to  numerically  determine  the 
domains  of  “growth  zone”  and  “stationary  zone”  for  a  given  tensile  stress  p  and  known  elastic 
material  properties.  Because  of  the  preferred  microcrack  growth  orientations,  the  overall  material 
compliances  become  even  more  anisotropic.  These  issues  will  be  further  addressed  in  Sec.  H.5. 

It  is  worth  noting  that  the  interfacial  fracture  criterion  for  a  tip  of  an  arc  microcrack  under 
biaxial  tension  field  is  also  available  from  Eq.  (4.7)  in  Toya  (1974).  Therefore,  the  regions  of 
microcrack  “growth  zone”  can  be  numerically  identified,  too. 
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As  previously  mentioned,  it  is  presumed  that  all  arc  microcracks  are  entirely  open  and  no  mi¬ 
crocrack  kinking  (into  matrix)  is  considered.  These  assumptions  are  made  because  solutions  for 
mixed  boundary-value  problems  are  not  yet  available  for  partially  closed,  entirely  closed  (mode 
II  frictional  sliding),  or  kinked  arc  microcracks.  When  these  solutions  do  become  available,  the 
proposed  process  model  must  be  properly  modified. 

In  addition,  it  is  recalled  that  the  region  of  “stable  microcrack  growth”  (from  an  initial  microc¬ 
rack  orientation  and  size  ( <f>it  a,)  to  a  final  stable  configuration  a })  is  quite  limited  and  is  only 

applicable  at  certain  critical  stress  level,  according  to  Fig.  3  in  Toya  (1974).  Therefore,  it  appears 
rational  and  practical  to  simplify  the  distributed  microcrack  growth  kinetics  as  follows.  If  one  tip 
of  an  arc  microcrack  reaches  or  exceeds  the  critical  surface  energy  required  to  initiate  microcrack 
growth,  then  both  tips  of  the  arc  microcrack  may  grow  continuously  (generally  in  a  non-symmetric 
fashion)  until  the  half-angle  size  reaches  an  a}  and  the  central  crack  ’’orientation”  tf>j  «  0°  (i.e., 
approximately  aligned  with  the  applied  tensile  loading  direction).  In  the  case  of  the  aforemen¬ 
tioned  epoxy-glass  composites,  o/  =  65°  is  a  reasonable  choice  since  it  is  the  maximum  allowable 
half-angle  size  for  an  arc  microcrack  to  remain  open .  It  is  noted  that  the  arc-microcrack  central 
orientation  indeed  will  change  due  to  crack  growth;  see  Toya  (1974)  for  more  discussions.  There¬ 
fore,  according  to  the  proposed  simplified  treatment,  an  arc  microcrack  may  either  be  stable  (no 
growth,  with  original  geometry)  or  unstable  (w,th  aj  =  65°  and  <f>j  =  0°).  Clearly,  damage-induced 
anisotropy  is  bound  to  be  a  natural  consequence. 

In  the  case  of  randomly  distributed  arc  microcracks,  Eq.  (4.2)  should  be  used  to  numerically 
(iteratively)  define  the  bounds  of  ( <f> ,  a)  regions  undergoing  microcrack  growth  under  a  particu¬ 
lar  tensile  stress  p,  a  given  interface  toughness  and  a  given  initial  microcrack  geometry.  For 
simplicity,  in  what  follows,  we  assume  that  all  arc  microcracks  are  of  equal  initial  size  a,  and  are 
such  oriented  that  \<t>\  +  a,  <  65°  (opening)  always  holds.  Thus,  for  a  specified  interface  toughness 
7*7,  one  can  perform  numerical  iteration  to  obtain  the  minimum  tension  stress  pCT  required  to  cause 
the  first  arc  microcracks  to  propagate.  The  corresponding  initial  central  microcrack  orientations 
are  denoted  by  ±<j>cr  .  When  p  is  greater  than  p„  ,  the  unstable  microcrack  domains  will  increase 
accordingly. 

The  proposed  microcrack  growth  kinetic  sequence  proceeds  as  follows. 
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(1)  As  p  <  per  ,  all  arc  microcracks  are  stable  and  of  initial  size  a,  .  Due  to  preferred  initial 
microcrack  orientations  to  ensure  opening,  the  overall  stress-strain  response  is  anisotropic  even 
if  virgin  composites  are  isotropic  (in  plane  strain).  Although  the  overall  response  is  linear  and 
reversible  under  the  present  stress  level  p ,  the  material  state  is  really  elastic-damage  (inelastic). 
In  fact,  the  current  elastic-damage  compliance  S  is  greater  than  the  virgin  undamaged  elastic 
compliance  S° . 

(2)  As  p  =  per ,  those  microcracks  with  central  orientations  <$>„  become  unstable,  and  increase  their 
lengths  to  aj  =  65°  as  well  as  change  their  orientations  to  <f>j  =  0 .  The  contributions  from 
partially  closed  and  entirely  closed  portions  are  neglected  here. 

(3)  As  p  =  pi  >  per  ,  the  unstable  microcrack  growth  regions  Qact  increase.  Therefore,  some 
microcracks  in  specific  orientation  domains  (-^/  ,  -<f>h)  and  ,  <f>i)  become  activated  and 
increase  in  size  to  ay  =  65°  as  well  as  change  their  orientations  to  <j>j  =  0  .  The  material 
behaves  anisotropically  and  the  elastic -damage  compliance  increases.  It  is  noted  that  4>i  and 
<f>h  values  depend  on  p\  ,  pcr  ,  a,  and  7'^,  and  can  be  obtained  by  numerical  iteration.  The 
compliance  contributions  S*  and  S*  in  Eq.  (41)  can  be  computed  (integrated)  through  Eq.  (35) 
and  (39): 


s;  =  N  /  g(k)Ts^kH<f>,  a f)gw  PWi> 

J  nacl 

(43) 

s*  =  N  [  gWTS *<*>'(<*,  a,)gw  Pm i 

Jn,„ 

(44) 

where  P(4> )  is  the  assumed  probability  density  function  of  microcrack  orientation  and  f2tn  is 
the  stable  (no  growth)  domain. 

(4)  As  pCT  <  p  <  pi  ,  the  unloading  case  is  taking  place.  There  is  no  further  microcrack  growth 
because  the  apparent  “active  orientation  fans”  shrink.  Therefore,  S*  =  0  .  It  is  emphasized 
that  the  actual  “orientation  fans”  (featuring  qj  size)  does  not  reduce  owing  to  the  irreversible 
nature  of  damage.  Hence,  the  elastic-damage  compliance  remains  its  previous  value. 

(5)  As  p  >  pi ,  more  microcracks  are  activated.  The  “unstable  growth”  domains  can  be  computed 
similar  to  step  (3).  However,  S*  in  step  (3)  should  be  replace  by  the  sum  of  S*  and  S}  in  this 
step. 

(6)  At  some  higher  stress  level  p  =  pc ,  the  energy  release  rate  reaches  the  critical  value  of  the 
matrix  energy  barrier.  Therefore,  microcracks  having  size  af  will  resume  to  propagate  (kink) 
into  the  matrix,  and  eventually  lead  to  final  failure. 
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As  was  commented  by  Sumarac  and  Krajcinovic  (1987),  the  above  scheme  implicitly  assumes 
that  ultimate  failure  prefers  “runaway  cracks”  in  comparison  with  “localization  modes”.  Some 
numerical  simulations  will  be  given  in  Sec.  II.5. 
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II.5.  Computational  algorithms  and  numerical  simulations 

In  this  section,  computational  integration  algorithms  are  given  for  the  proposed  micromechani¬ 
cal  damage  model.  Furthermore,  a  number  of  uniaxial  tension  numerical  simulations  are  presented. 
In  the  absence  of  suitable  plane  strain  experimental  data  at  this  stage,  however,  actual  experimental 
validation  is  not  included  here.  Experimental  verification  of  the  proposed  models  should  be  per¬ 
formed  in  the  future  once  they  become  available.  Nonetheless,  the  presented  numerical  simulations 
demonstrate  the  potential  capability  of  the  proposed  damage  model,  without  resorting  to  any  fitted 
“material  parameters”  commonly  utilized  in  phenomenological  continuum  damage  models. 

II.5.1.  Computational  integration  algorithms 

The  proposed  micromechanical  damage  model  does  not  include  microcrack  interaction  effects, 
and  therefore  falls  into  the  category  of  “Taylor’s  model”.  That  is,  one  employs  the  undamaged 
elastic  material  properties  (not  overall  elastic-damage  moduli)  to  compute  microcrack  opening  dis¬ 
placements.  Thus,  the  proposed  process  model  does  not  require  the  use  of  an  iterative  scheme  when 
computing  the  single  microcrack-induced  inelastic  compliance  S*(fc)  (see  Eqs.  (35H38)).  In  order 
to  determine  the  bounds  of  unstable  crack  growth  regions  under  a  particular  applied  tension  stress, 
nevertheless,  one  needs  to  perform  numerical  iterations,  as  commented  in  Sec.  II.4.2.  The  com¬ 
putational  schemes  involved  in  solving  the  proposed  stress-controlled  micromechanical  damage 
models  proceed  as  follows.  It  is  assumed  that  all  initial  arc  microcracks  are  of  equal  size  a,  and 
open. 

(1)  For  a  given  load  level  p  ,  compute  “unstable  orientation  bounds”  (- <f>t  ,  -<f>h)  and  (<^  ,  4>i) 
according  to  the  fracture  criterion  Eq.  (42).  For  example,  one  may  use  the  bi-section  method 
to  locate  the  very  first  unstable  microcrack  orientations  ±4>„  ,  and  later  the  bounds  defining 
the  crack  growth  domains.  These  “unstable  orientation  bounds”  should  be  stored  as  history 
variables  since  they  are  irreversible. 

(2)  Obtain  the  individual  and  total  damage-induced  inelastic  compliance  S*(i)  and  S*  by  actually 
evaluating  the  double  integral  involved  in  Eq.  (39)  (in  conjunction  with  Eqs.  (35)— (38)).  Note 
that  stable  microcTacks  have  the  size  a,  while  unstable  microcracks  reach  the  size  qj  (=  65° 
in  our  case).  Numerical  integration  of  Eq.  (39)  can  be  efficiently  performed  by  two  separate 
“Gauss  quadratures”  -  one  for  stable  and  the  other  for  unstable  regions.  In  particular,  one  needs 
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to  integrate  the  compliance  contributions  from  every  6  angle-coordinate  along  an  arc,  and  from 
every  open  arc  microcrack,  at  different  orientations.  Here,  we  use  Gauss  quadrature  rule  with 
60  integration  points  for  each  stable  or  unstable  region. 

(3)  Obtain  new  overall  elastic-damage  compliance  S  by  adding  S’  to  S° . 

(4)  Apply  the  next  load  p  and  go  to  step  (1). 

II.5.2.  Some  uniaxial  tension  tests 

A  number  of  mixed-mode  plane  strain  uniaxial  tension  tests  are  considered  in  this  section  for 
various  different  fiber  sizes  a,  microcrack  concentrations  u>  (=  N(aa)2/A),  and  initial  microcrack 
half-angle  sizes  a.  Both  stationary  and  process  (kinetic)  damage  models  are  employed  to  predict 
the  continuous  changes  in  overall  elastic-damage  compliances  and  stress-strain  responses  based  on 
micromechanics.  The  composite  material  considered  here  is  epoxy  for  matrix  and  glass  for  fiber. 
Both  epoxy  and  fibers  are  isotropic  and  linear  elastic.  The  shear  modulus  and  Poisson’s  ratio  of 
the  epoxy  matrix  are  p\  =  346  KSI  (2.39  GN/m2)  and  r/i  =  0.35,  respectively.  In  addition,  the 
shear  modulus  and  Poisson’s  ratio  of  the  glass  fibers  are  P2  -  6410  KSI  (44.2  GN/m2)  and  %  = 
0.22,  respectively.  The  composite  volume  ratio  of  the  matrix  and  fiber  constituents  are  0.8  and  0.2, 
respectively.  The  (plane  strain)  overall  composite  Young’s  modulus  and  Poisson’s  ratio  is  found  to 
be  E  =  1151  KSI  (7.93  GN/m2)  and  rj  =  0.3481,  respectively,  by  means  of  the  elementary  rule  of 
mixture  (see,  e.g.,  Jones  (1975)).  Therefore,  the  plane  strain  elastic  compliance  matrix  takes  the 
form:  (unit  =  KSI-1) 

‘  0.00076383  -0.00040784  0 

[S°]  =  -0.00040784  0.00076383  0  (45) 

0  0  0.0023433. 

Moreover,  the  interfacial  specific  surface  tension  energy  (fracture  toughness)  is  taken  as  0.001 
Kips/in.  (0.175  KN/m). 

First,  we  examine  the  effects  of  different  initial  microcrack  concentrations  uj,  on  stress-strain 
responses  and  compliances  of  the  uniaxially  reinforced  epoxy-glass  composite.  The  initial  micro¬ 
crack  concentration  parameter  u \  is  increased  gradually  (with  256  increments)  from  0  to  0.64.  All 
microcracks  are  assumed  to  be  stationary  (no  growth),  open,  and  of  the  half-angle  size  a,  =  10°. 
In  view  of  Fig.  2  in  Toya  (1974),  \<j>\  +  a  <  65°  is  required  for  crack  opening.  Therefore,  arc 
microcracks  are  assumed  to  be  perfectly  (uniformly)  randomly  oriented  between  4>  =  —55°  and  <f> 
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=  55°;  i.e.,  the  probability  density  of  orientation  is  P(4>)  =  0.521.  The  uniaxial  tension  stress  is 
applied  in  the  2-direction,  and  the  lateral  direction  is  denoted  as  the  1 -direction.  The  axial-stress 
(KSI)  versus  the  axial-  and  lateral-strain  (curves  corresponding  to  monotonically  increasing  values 
of  u>i  are  displayed  in  Fig.  2.  The  elastic  stress-strain  response  is  also  shown  in  Fig.  2  for  com¬ 
parison  purpose.  The  elasdc-damage  compliance  components  522  and  5i2  vs.  u>,  parameter  are 
shown  in  Fig.  3  and  Fig.  4,  respectively.  It  is  clear  that  as  microcrack  concentration  increases,  the 
compliance  components  522  and  5,2  increase  monotonically. 

Next,  we  perform  “process  model”  uniaxial  tension  tests.  Let  the  fiber  radius  be  a  =  0.1  in. 
and  microcrack  number  density  be  N/  A  =  100.  Two  initial  half-angle  sizes  of  arc  microcracks  are 
considered:  o,  =  10°  and  a,  =  20°,  respectively.  The  axial-stress  vs.  the  axial-  and  lateral-strain 
responses  are  recorded  in  Fig.  5.  It  is  noted  that  before  the  load  p  reaches  a  certain  critical  value  , 
the  stress-strain  response  is  linear  (stationary  damage  model)  for  either  a,  =  10°  or  a,  =  20°.  After 
the  critical  values  p„  are  reached,  the  overall  responses  then  become  nonlinear.  Double  Gauss 
quadrature  is  employed  to  integrate  inelastic  compliances  and  strains  contributions  from  various 
microcrack  orientations  and  sizes.  For  example,  in  the  case  of  a,  =  10°,  it  is  found  (by  numerical 
iteration)  that  the  first  microcracks  to  become  unstable  (propagating)  are  oriented  in  the  direction 
<t>c r  =  5.01°.  Later,  as  p  increases,  more  “orientation  fans”  are  enclosed  within  the  angle  bounds 
(obtained  by  numerical  iteration)  of  microcrack  growth.  Therefore,  more  and  more  microcracks 
experience  growth  and  finally  have  <f>j  -  0°  and  qj  =  65° .  In  the  case  of  a,  =  20°,  on  the  other 
hand,  it  is  found  that  (}>„  -  10.59°.  It  is  also  observed  from  Fig.  5  that  larger  initial  microcrack  size 
a,  =  20°  results  in  lower  critical  load  level  p„  for  microcrack  growth. 

Let  us  now  fix  the  initial  microcrack  size  a*  =  10°  and  vary  the  fiber  sizes  and  microcrack 
number  densities:  (i)  a  -  0.05  in.  anti  N/A  =  40,  and  (ii)  a  =  0.1  in.  and  N/A=  100,  respectively. 
The  axial-stress  vs.  the  axial-  and  lateral-strain  responses  are  shown  in  Fig.  6.  The  changes  in 
overall  elastic -damage  compliances  522,  •S’,,  and  5,2  versus  axial  stresses  are  plotted  in  Fig.  7,  8 
and  9.  It  is  observed  that  522  and  5,,  increase  as  tensile  stress  level  p  increases.  However,  5,2  is 
negative  and  decreases  as  the  stress  increases;  i.e.,  |5,2|  increases  monotonically. 

The  above  uniaxial  tension  tests  are  repeated  for  fixed  initial  microcrack  size  a,  =  20°.  Again, 
two  sets  of  fiber  sizes  and  microcrack  number  densities  are  considered:  (i)  a  =  0.05  in.  and  NjA  = 
40,  and  (ii)  a  =  0.1  in.  and  N/A  =  100.  The  axial-stress  vs.  the  axial-  and  lateral- strain  responses 
are  shown  in  Fig.  10.  The  changes  in  overall  elastic-damage  compliances  S22,  5,,  and  5,2  versus 
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axial  stresses  are  plotted  in  Fig.  1 1, 12  and  13.  Similarly,  S&  and  Sn  increase  as  tensile  stress  level 
p  increases.  In  addition,  S n  decreases  as  the  stress  increases. 

The  above  numerical  simulations  illustrate  the  potential  capability  of  the  proposed  microme¬ 
chanical  damage  models  to  qualitatively  explain  and  quantitatively  predict  the  overall  stress-strain 
behavior  and  compliance  evolutions.  Not  a  single  fitted  “material  parameter”  is  employed  in  the 
proposed  model.  It  is  also  emphasized  that  the  numerical  tests  presented  in  this  section  is  for 
“cleavage  1”  damage  processes.  Thus,  nucleation  of  new  arc  microcracks  is  not  included.  Fur¬ 
ther,  partially  or  entirely  closed  arc  microcracks,  or  microcrack  kinking  (into  the  matrix)  are  not 
incorporated  into  the  present  model.  These  issues  will  be  topics  of  future  research. 
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II.6.  Conclusion 

A  micromechanical  damage  model  is  presented  in  this  work,  based  on  Toya’s  (1974)  microme¬ 
chanics  solutions  and  general  damage  mechanics  framework,  for  uniaxially  reinforced  brittle  ma¬ 
trix  fibrous  composites.  All  microcracks  are  assumed  to  exist  along  the  fiber/matrix  interfaces,  and 
are  modeled  as  arc  microcracks  under  “cleavage  1”  plane  strain  deformation  processes.  Thermody¬ 
namic  basis  is  rendered  based  on  a  fourth-rank  damage  tensor.  Microcrack-induced  inelastic  strains 
and  elastic-damage  compliances  are  analytically  derived  based  on  micromechanical  bi-material  arc- 
microcrack  opening  displacements  and  probabilistic  distributions  of  arc  microcracks.  It  is  noted 
that  the  (plane  strain)  overall  elastic-damage  compliance  matrix  is  non-symmetric.  Both  “ station¬ 
ary ”  and  “ process ”  damage  models  are  given  in  the  present  work.  In  particular,  microcrack  kinetic 
equations  are  given  based  on  micromechanical  fracture  criterion  and  microstructural  geometry  in 
a  representative  volume  (area)  element.  Moreover,  simple  computational  algorithms  and  a  number 
of  uniaxial  tension  tests  are  presented  to  illustrate  the  potential  capability  of  the  proposed  microme¬ 
chanical  damage  model  for  fibrous  composites.  It  is  emphasized  that  the  present  work  does  not  em¬ 
ploy  any  arbitrary  (fitted)  “material  constants”.  The  proposed  framework  can  be  readily  extended 
to  account  for  biaxial  tension  loadings,  as  addressed  in  Remark  3.1.  Issues  related  to  microc¬ 
rack  interaction,  closed  microcracks,  microcrack  kinking,  and  microcrack  nucleation  mechanisms 
warrant  further  studies  in  the  future  to  extend  the  proposed  method. 
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II.8.  Figure  captions 

Figure  1.  The  local  (primed)  and  global  Cartesian  coordinates,  as  well  as  the  polar  coordinates 
at  a  typical  arc  point. 

Figure  2.  The  axial-stress  vs.  the  axial-  and  lateral- strain  under  different  microcrack  concen¬ 
trations  u  (varying  from  0  to  0.64).  Note  that  u  =  N(aa)2/ A .  The  dotted  lines  are  the  undamaged 
elastic  stress-strain  responses. 

Figure  3.  The  overall  compliance  S a  vs.  microcrack  concentration  u>  for  various  stationary 
damage  model  simulations.  The  dotted  line  is  the  elastic  response. 

Figure  4.  The  overall  compliance  Sn  vs.  w  for  various  stationary  damage  model  simulations. 

Figure  5.  The  axial-stress  vs.  the  axial-  and  lateral-strain  for  two  different  initial  microcrack 
sizes:  a,  =  10°  and  a,  =  20°.  Note  that  a  =  0.1  in.  and  N/A  =  100.  These  are  “process  model” 
simulations. 

Figure  6.  The  axial-stress  vs.  the  axial-  and  lateral-strain  for  two  different  fiber  sizes  and 
microcrack  number  densities.  Note  that  a,  =  10°. 

Figure  7.  The  overall  compliance  S22  vs.  the  axial  stress  p  for  two  different  sets  of  a  and  N/A 
values.  Note  that  a;  =  10°. 

Figure  8.  The  overall  compliance  Sn  vs.  the  axial  stress  p  for  two  different  sets  of  a  and  N/A 
values.  Note  that  Oi  =  10°. 

Figure  9.  The  overall  compliance  vs.  the  axial  stress  p  for  two  different  sets  of  a  and  N/A 
values.  Note  that  a,  =  10°. 

Figure  10.  The  axial-stress  vs.  the  axial-  and  lateral- strain  for  two  different  fiber  sizes  and 
microcrack  number  densities.  Note  that  a,  =  20°. 

Figure  11.  The  overall  compliance  S&  vs.  the  axial  stress  p  for  two  different  sets  of  a  and 
N/A  values.  Note  that  a,  =  20°. 

Figure  12.  The  overall  compliance  Sn  vs.  the  axial  stress  p  for  two  different  sets  of  a  and 
N/A  values.  Note  that  a,  =  20°. 

Figure  13.  The  overall  compliance  S12  vs.  the  axial  stress  p  for  two  different  sets  of  a  and 
N/A  values.  Note  that  a,  =  20°. 


Figure  2.  The  axial-stress  vs.  the  axial-  and  lateral-strain  under  different  microcrack  concen¬ 
trations  u>  (varying  from  0  to  0.64).  Note  that  u  =  N(aa)2/A .  The  dotted  lines  are  the  undamaged 
elastic  stress-strain  responses. 
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Figure  3.  The  overall  compliance  S22  vs.  microcrack  concentration  u  for  various  stationary 
damage  model  simulations.  The  dotted  line  is  the  elastic  response. 


(i-ISX)  £-01  X80irei[dui03 


DAMAGED 

ELASTIC 


Figure  4.  The  overall  compliance  Sn  vs.  u  for  various  stationary  damage  model  simulations. 
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Figure  5.  The  axial-stress  vs.  the  axial-  and  lateral-strain  for  two  different  initial  microcrack 
sizes:  a,  =  10°  and  a,  =  20°.  Note  that  a  =  0.1  in.  and  N/A  =  100.  These  are  “process  model” 
simulations. 
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Figure  6.  The  axial-stress  vs.  the  axial-  and  lateral-strain  for  two  different  fiber  sizes  and 
microcrack  number  densities.  Note  that  a;  =  10°. 
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Figure  7.  The  overall  compliance  S2 2  vs.  the  axial  stress  p  for  two  different  sets  of  a  and  N/A 
values.  Note  that  a,  =  10°. 
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Figure  8.  The  overall  compliance  Sn  vs.  the  axial  stress  p  for  two  different  sets  of  a  and  N/A 
values.  Note  that  a,-  =  10°. 
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Figure  9.  The  overall  compliance  S12  vs.  the  axial  stress  p  for  two  different  sets  of  a  and  N/A 
values.  Note  that  a,-  =  10°. 
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Figure  11.  The  overall  compliance  S22  vs.  the  axial  stress  p  for  two  different  sets  of  a  and 
N/A  values.  Note  that  a,-  =  20°. 
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A  Three-Dimensional  Statistical  Micromechanical 


Theory  for  Brittle  Solids  with  Interacting  Microcracks 

III.O.  Abstract 

A  three-dimensional  statistical  micromechanical  theory  is  presented  to  investigate  the 
overall  nonlinear  mechanical  responses  of  microcrack-weakened  brittle  solids.  The  macro¬ 
scopic  stress-strain  relations  of  elastic  solids  with  interacting  microcracks  are  micromechan- 
ically  derived  by  taking  the  ensemble  average  over  all  possible  realizations  which  feature  the 
same  material  mesostructural  geometry.  Approximate  analytical  solutions  of  a  two-crack  in¬ 
teraction  model  are  introduced  to  account  for  microcrack  interaction  among  many  random¬ 
ly  oriented  and  distributed  microcracks.  The  overall  effective  compliances  of  microcrack- 
weakened  brittle  solids  are  derived  by  further  taking  the  volume  average  of  the  ensemble- 
averaged  stress-strain  relations  over  the  entire  material  mesostructural  domain  of  a  represen¬ 
tative  volume  element.  Moreo\  er,  some  special  cases  are  investigated  by  using  the  proposed 
framework  and  other  existing  methods.  It  is  emphasized  that  no  Monte  Carlo  simulations 
are  needed  in  the  proposed  framework. 
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The  nonlinear  mechanical  responses  of  damaged  solids  due  to  the  existence,  growth,  and 
nucleation  of  microdefects  (such  as  microcracks  and  microvoids)  are  of  significant  importance 
to  engineers,  and  have  been  the  subject  of  many  investigations.  See  Krajcinovic  (1989)  for 
an  excellent  literature  review  on  damage  mechanics.  For  brittle  materials  (e.g.  concrete, 
rocks  and  ceramics),  in  particular,  microcracks  often  control  overall  deformation  and  failure 
mechanisms.  To  date,  the  only  exact  results  derived  for  microcrack-weakened  brittle  solids 
are  for  dilute  microcrack  concentrations,  where  microcrack  interactions  are  entirely  neglected. 
These  micromechanical  damage  models  are  called  “Taylor’s  models”;  see,  e.g.,  Krajcinovic 
and  Fanella  (1986),  Fanella  and  Krajcinovic  (1988),  and  Ju  (1990)  for  some  recent  works. 

On  the  other  hand,  several  approximate  micromechanical  analyses  (“effective  medium 
methods”)  were  proposed  in  the  literature  to  account  for  interaction  effects  of  distributed 
microcracks.  For  example,  the  “self-consistent  method”  (Hill  (1965))  was  first  applied  to 
microcrack-weakened  solids  by  Budiansky  and  O’Connell  (1976)  with  special  attention  di¬ 
rected  to  perfectly  randomly  distributed  (isotropic)  and  weakly  interacting  microcracks.  The 
self-consistent  method  was  further  developed  by  Horii  and  Nemat-Nasser  (1983)  to  take  into 
account  the  effects  of  closed  microcracks  undergoing  frictional  sliding.  See  also  Sumarac  and 
Krajcinovic  (1987,  1989),  Krajcinovic  and  Sumarac  (1989),  Ju  (1991),  Ju  and  Lee  (1991), 
and  Lee  and  Ju  (1991).  Christensen  and  Lo  (1979)  proposed  a  three-phase  “genei  alized 
self-consistent  model”.  The  “differential  scheme”  was  investigated  by  Roscoe  (1952,  1973), 
McLaughlin  (1977),  and  Hashin  (1988).  Further,  the  “Mori-Tanaka  method”  was  developed 
by  Mori  and  Tanaka  (1973),  Benveniste  (1986),  and  Zhao,  Tandon  and  Weng  (1989).  Based 
on  variational  principles,  Hashin  and  Shtrikman  (1962,  1963)  also  proposed  upper  and  lower 
bounds  for  composites  with  inclusions. 

Some  comparisons  and  assessments  for  the  self-consistent  method,  the  generalized  self- 
consistent  method,  the  Mori-Tanaka  method,  and/or  the  differential  scheme  were  also  pre¬ 
sented  by  Horii  and  Sahasakmontri  (1990),  Laws  and  Dvorak  (1987),  Nemat-Nasser  and  Hori 
(1990),  and  Christensen  (1990,  for  pure  shear  load  only).  It  is  noted  that  effective  medium 
methods  are  only  valid  for  weak  or  at  most  moderate  microcrack  concentrations.  Overall 
effective  compliances  as  well  as  anisotropy  due  to  microcracks  can  be  estimated  by  effective 
medium  approximations. 
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When  microcrack  concentrations  are  higher  and  microcrack  spacings  are  closer,  strong 
microcrack  interactions  occur  and  effective  medium  theories  are  no  longer  appropriate.  Em¬ 
anating  from  this  viewpoint,  excellent  strong  microcrack  interaction  models  were  proposed 
by  Horii  and  Nemat-Nasser  (1985),  Hori  and  Nemat-Nasser  (1987)  for  two-dimensional  de¬ 
terministic  microcracks  (not  at  the  overall  constitutive  level),  and  by  Kachanov  and  Mon- 
tagut  (1986),  Kachanov  (1987),  Chudnovsky  et  al.  (1987a, b),  Montagut  and  Kachanov 
(1988),  Kachanov  and  Laures  (1989),  and  Laures  and  Kachanov  (1990)  for  two-  and  three- 
dimensional  deterministic  (specified)  arbitrary  microcrack  arrays  (at  the  overall  constitutive 
level).  Other  literatures  are  cited  and  discussed  in  the  reference  papers.  Moreover,  the 
important  and  valuable  work  due  to  Kachanov  (1987),  Kachanov  and  Laures  (1989)  rely 
on  Monte  Carlo  simulations  of  a  large  amount  of  deterministic  microcrack  arrays,  and  de¬ 
pend  heavily  on  large-scale  computations  to  obtain  local  stresses  of  many  randomly  located 
microcracks.  Therefore,  it  is  desirable  to  develop  innovative  and  simple  statistical  microme¬ 
chanical  damage  theories  to  account  for  interactions  among  randomly  located  and  oriented 
microcracks  without  using  extensive  random  microcrack  simulations  and  large-scale  iterative 
stress  computations. 

The  purpose  of  the  present  work  is  to  establish  a  three-dimensional  statistical  microme¬ 
chanical  framework  to  predict  overall  effective  moduli  for  brittle  solids  with  many  interacting, 
randomly  distributed  microcracks.  The  proposed  statistical  framework  considers  the  proba¬ 
bility  and  conditional  probability  density  functions  of  microcrack  locations  and  relative  con¬ 
figurations.  See  Batchelor  (1970),  Hinch  (1977),  Willis  (1977),  Chen  and  Acrivos  (1978a, b) 
for  references.  In  addition,  the  ensemble  and  volume  averages  of  stresses,  strains  and  com¬ 
pliances  are  systematically  constructed  based  on  analytical  micromechanics  solutions  and 
statistical  mechanics  concepts.  Therefore,  statistical  aspects  are  naturally  embedded  into 
damage  theories  and  macroscopic  (averaged)  constitutive  equations.  As  a  result,  we  do  not 
rely  on  heavy  numerical  stress  computations  or  extensive  Monte  Carlo  random  simulations 
involving  hundreds  of  deterministic  microcrack  arrays. 

It  is  emphasized  that  the  proposed  method  is  fundamentally  different  from  that  pro¬ 
posed  by  Hudson  (1980,  1981,  1986).  Hudson’s  method,  though  using  the  ensemble  average 
approach,  is  based  on  a  second  order  stiffness  theory  and  thus  leads  to  irrational  behaviors 
for  microcracked  solids  with  moderate  or  high  microcrack  concentrations  (see  Sayers  and 
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Kachanov,  1991).  That  is,  Hudson’s  method  is  designed  for  small  microcrack  density.  The 
proposed  approach  is  suitable  for  low,  moderate  or  high  microcrack  concentrations. 

A  brief  outline  of  this  work  is  as  follows.  In  section  2,  an  ensemble  average  approach  to 
derive  damaged  stress-strain  relations  is  introduced.  Approximate  closed-form  analytical  so¬ 
lutions  are  subsequently  presented  for  the  interaction  problem  of  two  aligned  but  arbitrarily 
located  penny-shaped  microcracks.  The  overall  moduli  of  a  representative  volume  element 
(unit  cell)  are  then  derived  by  further  taking  the  volume  average  of  ensemble  averaged  local 
stress-strain  relations  over  the  entire  (unit  cell)  domain  at  the  mesostructural  level  in  Sec¬ 
tion  3.  In  Section  4,  applications  are  made  to  a  number  of  special  cases.  In  particular,  for 
the  dilute  microcrack  concentration  case,  the  present  approach  recovers  the  well-known  Tay¬ 
lor’s  model  by  neglecting  interactions  among  microcracks.  Further,  for  aligned  penny-shaped 
microcracks,  the  present  method  reveals  that  average  tractions  and  overall  compliances  are 
changed  by  the  presence  and  interaction  of  3-D  microcracks.  This  result  is  relevant  to  fiber 
breaking  in  unidirectionally  reinforced  fiber  composites.  Comparisons  with  the  Taylor’s 
model,  the  self-consistent  method  and  the  differential  scheme  are  also  presented.  Finally, 
higher-order  microcrack  interaction  models  within  the  proposed  framework  are  discussed  in 
Section  5. 
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III.2.  An  ensemble  average  approach  to  3-D  aligned  microcrack 
interaction  and  effective  moduli 

We  first  review  some  basic  background  related  to  the  concept  of  volume-averaged  (mesostruc- 
tural)  stress  and  strain  tensors  within  the  context  of  a  representative  volume  element  (RVE) 
or  unit  cell.  The  ensemble  average  of  microcrack-induced  strains  and  approximate  analytical 
solutions  of  two  arbitrarily  located  and  oriented  microcracks  are  subsequently  given.  Overall 
effective  compliances  of  brittle  solids  with  interacting  microcracks  are  then  derived  within 
the  framework  of  statistical  micromechanics.  Following  current  literatures,  it  is  assumed 
that  the  volume- average  stress  &  approximately  equal  to  the  far- field  stress  <r°°.  This  typ¬ 
ical  assumption  may  be  removed  and  significant  further  improvements  are  very  promising. 
These  issues  will  be  addressed  in  a  forthcoming  paper. 

III.2.1.  Ensemble  average  of  microcrack-perturbed  stresses  and  s- 
trains 

Due  to  the  existence  and  interaction  of  microcracks,  local  stresses  and  strains  in  the 
matrix  material  are  perturbed.  The  ensemble  average  approach  hinges  on  the  concept  that 
local  stresses,  strains  and  compliances  (or  stiffnesses)  at  a  typical  point  within  a  RVE  of  a 
microcrack-weakened  solid  can  be  obtained  by  averaging  over  the  ensemble  of  all  statistical 
realizations  of  randomly  distributed  microcracks.  Batchelor  (1970),  Batchelor  and  Green 
(1972),  and  Hinch  (1977)  applied  this  approach  to  the  study  of  fluid  suspensions  within 
the  framework  of  pairwise  (second  order)  interaction.  The  ensemble  average  approach  with 
pairwise  interaction  w’as  later  applied  to  composite  materials  with  interacting  inclusions  (in¬ 
homogeneities)  by  Willis  and  Acton  (1976),  as  well  as  Chen  and  Acrivos  (1978a, b).  Recently, 

Ju  and  Chen  ( 1990a, b)  proposed  two-dimensional  (second  order  and  higher)  micromechan¬ 
ical  damage  theories  for  brittle  solids  with  interacting  slit  microcracks  by  employing  the 
ensemble  average  method  and  micromechanical  fracture  mechanics.  It  is  noted  that  local 
displacements,  strains  and  stresses  vary  with  positions  within  a  RVE. 

For  simplicity,  we  consider  a  two-phase  composite  composed  of  a  linear  elastic  matrix 
and  many  penny-shaped  microcracks.  The  local  strain  tensor  at  a  point  x  within  the  RVE 
takes  the  form: 


e(x)  =  S°  :  <r(x)  +  e*(x,C) 


(1) 
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where  e  and  a  denote  the  local  strain  and  stress,  respectively;  e*  is  the  perturbed  strain 
due  to  the  existence  and  interaction  of  microcracks;  and  C  denotes  the  set  of  all  possible 
configurations  of  the  microcracks.  It  is  emphasized  that  em  is  zero  if  x  is  a  point  in  the 
matrix  and  non-zero  if  x  is  a  point  on  the  microcrack  surfaces.  Taking  the  ensemble  average 
over  Eq.  (1),  we  arrive  at 

<«)(*)  =  S*  :  (<r)(x)  +  <OM  (2) 

where  the  angle  brackets  (•)  signify  the  ensemble  average. 

Throughout  the  development  of  this  paper,  for  simplicity,  we  shall  assume  that  the  solid 
is  locally  homogeneous  (Hinch  (1977))  and  penny-shaped  microcracks  do  not  intersect  one 
another.  Locally  homogeneity  implies  that  all  probability  density  functions  (PDF)  do  not 
vary  under  small  translation  on  a  macroscopic  length  scale.  It  can  be  shown  (Ju  and  Chen 
(1990a))  that,  in  the  case  of  microcracks,  the  perturbed  strain  can  be  expressed  as 

<*•)(*)  =  /(*)jfjf  j(H  ®  n  +  n  ®  |u])(x'|x,5)/(e)  dQiS  (3) 

Here,  Q  =  (a,  n)  indicates  the  microcrack  length  (radius)  a  and  orientation  n.  Furthermore, 
x'  denotes  a  point  on  the  surfaces  (5.)  of  a  microcrack  centered  at  x;  /(x)  is  the  PDF  for  a 
microcrack  being  centered  at  x;  Ju]  is  the  vector  of  microcrack  opening  displacements;  and 
f{Q)  is  the  PDF  for  a  microcrack  with  a  geometry  Q.  If  all  penny-shaped  microcracks  are 
aligned  (parallel)  and  of  equal  size,  then  there  is  no  variation  in  Q.  This  could  correspond  to, 
for  example,  microcracks  generated  by  fiber  breaks  in  unidirectional  fiber  composites  (see, 
e.g.,  Laws  and  Dvorak  (1987)).  In  this  event,  Eq.  (3)  can  be  simplified  as  follows 

(O(x)  =  /(x)  f  5«u]  ®  n  +  n  ®  Ju])(x/|x)  dS  (4) 

J  5, 


It  is  well  known  that,  for  an  open  penny-shaped  microcrack  with  radius  a  embedded  in 
an  infinite  linear  elastic  isotropic  matrix,  the  microcrack  opening  displacements  at  x'  (at  a 
distance  p  from  the  center  of  the  microcrack)  are: 


8(1  -  »7) 

irE(  2  —  v ) 


(5) 


where  E  and  v  =  the  Young’s  modulus  and  Poisson’s  ratio  of  the  virgin  matrix  material, 
respectively.  Moreover,  p,  s  and  t  =  the  z-direction  normal,  the  x-direction  shear  and  the 
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y-direction  shear  stresses  projected  on  the  microcrack  surface  in  its  local  coordinates;  see 
Fig.  1  for  a  schematic  plot. 


If  all  (open)  microcracks  are  aligned  (parallel),  then  we  can  define  n  =  (0, 0, 1)T.  By 
substituting  Eq.  (5)  into  (4)  and  earring  out  the  integration,  we  obtain  (assuming  equal 
microcrack  size) 

<*■><*>  =  8 ' (T>  (6) 

where  (the  Voigt’s  notation) 


(7) 


In  addition,  g  is  the  transformation  matrix  and  T  is  the  local  stress  vector: 


'0  0  0‘ 

0  0  0 

„  _  2-i/  0  0 

g~  0  0  0 

0  0  2 

0  2  0 


(8) 


In  the  event  of  distributed  (nonuniform)  microcrack  lengths  and  orientations,  similar  expres¬ 
sions  can  be  constructed  accordingly. 


On  the  other  hand,  the  local  stress  vector  (T)  can  be  shown  to  be  (see  Eq.  (14)  in 
Sec.  2.2): 

T  =  T00  +  T=  js00J  +  jfj 

Here,  T°°  denotes  the  unperturbed  local  stress  vector  due  to  remote  loading,  and  T  denotes 
the  local  stress  perturbation  due  to  three-dimensional  microcrack  interactions.  In  what 
follows,  attention  is  focused  on  pairwise  (second  order)  microcrack  interactions.  Higher 
order  microcrack  interactions  will  be  discussed  in  Sec.  5.  By  assuming  that  all  microcracks 
are  aligned  with  a  chosen  global  coordinate  system,  the  stress  T°°  due  to  far  field  loads  can 
be  expressed  as 


T°°  =  K0  t°° 


(10) 


A  3-D  Statistical  Micromechanical  Theory 


59 


where 


K0 


'0  0  1  0  0  0' 

0  0  0  0  0  1 

0  0  0  0  1  0 


XX 

T°° 

yy 

-oo 


zy 

r°° 

'yz 

Too 


(11) 


The  ensemble  average  of  local  stress  perturbation, 


on  the  other  hand,  takes  the  form  (as¬ 


suming  uniform  size  and  aligned  orientation): 


(T)(x)  =  j((T)(x|x,)/(x,|x)  dx, 


(12) 


where  (T)(x|xi)  =  the  ensemble- average  stress  perturbation  for  a  microcrack  centered  at 
x  over  the  subclass  of  realizations  having  a  microcrack  centered  at  Xj;  and  /(xj|x)  =  the 
conditional  PDF  for  finding  a  microcrack  centered  at  xx  given  a  microcrack  centered  at  x. 
Further,  E  designates  the  active  (open)  integration  domain  which  depends  on  the  loading 
conditions. 


The  conditional  PDF  /(x^x)  can  be  simplified  to  /(xi)  if  microcracks  do  not  intersect 
and  reasonable  randomness  holds  (Hinch  (1977)).  The  local  homogeneity  assumption  enables 
us  to  further  approximate  /(xj)  by  /(x).  Therefore,  Eq.  (12)  can  be  recast  as: 

(T)(x)  =  /(x)  /(t)(x|x,)  At,  (13) 

The  quantity  T(x(xi)  corresponding  to  pairwise  microcrack  interaction  will  be  the  main 
subject  of  Sec.  2.2. 


III.2.2.  Approximate  explicit  solutions  for  pairwise  interaction  of 
aligned  microcracks 

The  objective  of  this  section  is  to  construct  approximate  closed-form  explicit  expressions 
for  perturbed  stresses  T  due  to  two-microcrack  interaction  so  that  the  ensemble- average 
formalism  proposed  in  Sec.  2.1  can  be  realized.  It  is  possible,  in  principle,  to  derive  expres¬ 
sions  for  T  in  the  interaction  problem  involving  two  3-D,  arbitrarily  oriented  and  located 
penny-shaped  microcracks.  However,  the  general  solutions  of  two-microcrack  interaction  are 
rather  complicated  and  no  reasonably  compact  closed-form  explicit  solutions  are  possible 
(although  numerical  solutions  are  certainly  feasible).  Therefore,  for  demonstration  purpose, 
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attention  will  be  foucused  on  explicit  solutions  of  two  randomly  located,  aligned  (parallel), 
equally  sized,  penny-shaped  microcracks  embedded  in  infinite  linear  elastic  isotropic  matrix. 

The  “pseudo-traction”  method  is  adopted  here  to  derive  approximate  expressions  for  T 
since  exact  solutions  are  not  yet  available.  For  mathematical  simplicity,  only  the  first  term 
of  Taylor’s  expansion  of  the  local  stress  field  is  used  to  represent  the  average  stress  across  the 
microcrack  surface.  Higher  order  terms  in  polynomial  expansions  may  be  included  if  desired, 
however,  at  the  high  cost  of  a  much  larger  system  of  equations  and  much  more  complicated 
analytical  expressions.  Stemming  from  a  different  viewpoint,  the  more  accurate  “transmis¬ 
sion  factor”-type  formulation  proposed  by  Kachanov  (1987)  is  well  suited  for  deterministic, 
numerical  computations.  Within  the  framework  of  statistical  ensemble  average,  nevertheless, 
closed-form  explicit  solutions  for  stress-interaction  are  warranted. 

Figure  1  shows  the  local  coordinate  systems  for  microcracks  1  and  2  of  radius  a.  The 
2-axis  is  chosen  as  the  direction  normal  to  a  microcrack  surface.  In  accord  with  the  pseudo¬ 
traction  concept,  the  problem  of  two  interacting  microcracks  subjected  to  far  field  stresses 
can  be  decomposed  into  a  homogeneous  problem  and  two  sub- problems  (see  also  Ju  and 
Chen  (1990a)'l.  Tn  the  homogeneous  problem,  a  microcrack-free  solid  is  subjected  to  applied 
stresses  at  far  field.  In  the  sub-problem  j  ( j  =  1,  2),  an  infinitely  extended  solid  contains  only 
one  penny-shaped  microcrack  and  is  subjected  to  zero  remote  stress  at  infinity.  Since  stresses 
vanish  on  microcrack  surfaces  (Cj),  the  following  boundary  conditions  must  be  satisfied  ( j 
=  1,  2): 

~Pj  +  +  Pj  —  0  •  —Sj  +  s°°  +  Sj  =  0  ;  —tj  +  t  °°  +  ii  =  0  (14) 

In  the  first  sub-problem  (containing  only  ‘microcrack  1’),  let  us  define  a  cylindrical 
coordinate  system  with  the  center  of  the  ‘microcrack  1’  as  its  origin.  Accordingly,  the  center 
location  of  the  ‘microcrack  2’  can  be  characterized  by  (/?,  <f>,  2);  see  Fig.  1.  For  convenience, 
we  shall  introduce  the  following  definitions: 

O  \  ~  G XX  "1"  Gyy  ,  (T  2  ==  CT  XX  CFyy  2  2  (Jj-y  (15) 

(J z  —  O zz  j  rx  —  &xz  -f"  i&yz  (1®) 

If  the  surfaces  of  the  ‘microcrack  1’  are  subjected  to  applied  normal  stress  p^\  then  perturbed 


A  3-D  Statistical  Micromechanical  Theory 


61 


stresses  at  the  ‘microcrack  2’  location  are  given  by  (Fabrikant  (1989,  p.  252-257)): 


«?>  -  ((1 + 8,)  \a{%;  -  si„-  ( f ' )]  +  «,w+«,(a.,+fa,-3^)ii 

*  V  [  Q-<i  Wj  (<1  -  If)'  01  -  a2)  I 

m  _  2p'»  ajy^q  -  a2)'/2  r  _  ^  [»a(ei|  -  2 1;  +  Pp  -  sij)  1 

2  »  m-t)  \  1  w-mii-a’)  } 

41)  2p(l)  (  a(l(  -  a2y/'!  .  / a\  az2  [/)  +  a1  (‘la1  +  2z2  -  3/)2 )]  ] 

'  "  •  1  <?-<?  Sm  U/  (?  -  W(P  -  J 

4.)  _  2pW  [a;(4lj  -  5p»)  +  If) 

W'l  -  l})3 

where 

^  =  2  {  [(a  +  ^  +  ^2]1/2  ~  [(a  ~  ^)2  +  z2]1/2} 

^2  =  2  {  [(a  +  ^)2  +  ^2]1/2  +  [(a  -  p)2  +  z2]1/2| 


Furthermore,  when  the  surfaces  of  the  ‘microcrack  1’  are  loaded  by  the  combined  shear 
stresses  then  perturbed  stresses  at  the  ‘microcrack  2’  location  are  (Fabrikant 

(1989,  p.  257-261)): 

2(f<»e'*  +  T<»e-**)  f  *1,^ -«>)■/» [«»(4I|- */>>)  + /})! 

1 - ?(2^)  r2(1  +  i,)  wF/j)  +  5iRf  / 


»  _  -2^  f4H  _  a zt,(g  -  <■»)'/»  f4£  (1)  2rt  _  c*(4Jj  -  V)  + ;;  „J  1 

1  t(2  -  p)  \  11  '  Wl  -12)  +  /2(ij  -ID  [II  (II  -  ID 2  1  J  J 


2(f0>e«  +  rC>e-'D  zl,(l2  -  a2)2'2  [aJ(4/§  -  5p!)  +  IJ] 

1  »(2  -  p)  wi  -  i2,r 


_  2  /  L  _  -  °2)1/;  _  .  ,  Z{ °2  -  f?)l/;  W  +  42(2qj  +  2z;  -  3p2)] 

x(2-p)\r  'V  2|  -  ^  UJJ  01 -‘I)3 

,  L,,*  *(»2  -  ID'12  [«2(642  -  21*  +  p*)  -  5f<)l  IJ«W  ml 

+  rft-*)  + - oFifF - / 


where  r*1)  is  the  complex  conjugate  of  r^.  The  total  perturbed  stresses  due  to  combined 

normal  loading  pM  and  shear  loading  are  therefore  obtained: 

-0)  _  _/(»)  ,  _//(i)  .  _(i)  _  _/(i)  ,  _«(i) 

"■  -<'1  +<Tl  ’  +,T2  (20) 

<r<”  =  +  <r"i”  ;  r<>»  =  r'i11  +  r"'1’ 

It  can  be  shown  that  the  Cartesian  stress  components  are 

«£  =  5  d”  +  R*^11)]  ;  <*11’  =  5  p!11  -  Re^11)] 

*£?-• f>  i  ^  =  <21) 

=  Im(rz(1))  ;  o{£  =  Im(rz(,)) 
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where  “Re”  and  “Im”  are  the  real  and  imaginary  parts,  respectively,  of  a  complex  variable. 


Eq.  (21)  can  be  recapitulated  into  the  following  matrix  form: 

/  \  (1)  r  ti  +b7  c,+C2  dj+di  -1  (*) 


r'  '  = 


_  \ 
0 XX 

(1) 

Oyy 

<?zz 

>  ~ 

0Xy 

&yz 

2  2  2 
61—62  Cl  *— C2  rfl 


a 

,2  2 

t>5  C5 

be  ce 


(3' 


where  definitions  of  the  parameters  6,,  c^,  and  d<  can  be  found  in  Appendix  I.  Symbolically, 
we  can  express  as  a  function  of  coordinates  and  normal  and  shear  stresses  as  follows: 

t(1)  =  r(p,4>,2-,p1,s1,t1)  (23) 


Since  microcracks  1  and  2  are  aligned,  the  perturbed  normal  and  shear  stresses  along 
the  surfaces  at  the  ‘microcrack  2’  location  are 

P2  =  e2  •  r(1)  •  e2  = 

h  =  e*  •  r(1)  •  ex  =  (24) 

t2  =  e2  •  •  ev  =  e™ 

where  ex,  ey,  and  ez  are  the  unit  base  vectors  in  the  Cartesian  (x,  y,  and  z)  coordinates. 

Similarly,  in  the  sub-problem  2  (containing  only  ‘microcrack  2’),  it  can  be  shown  that 
the  perturbed  normal  and  shear  stresses  along  the  surfaces  at  the  ‘microcrack  1’  location  are 

Pi  =  vi?  ;  ^  =  -°ix  i  *1  =  0{$  (25) 

where  symbolic  representations  similar  to  Eq.  (23)  and  (24)  have  been  employed: 

t{2)  =  r{p,n  -  <f>,z-,p2,-S2,t2)  (26) 

Eq.  (24)  together  with  (25)  then  leads  to 

°  0  0  6<J)  -4!)  d?1  1 

?;  000  -if  4”  -42)  5] 

|  <■  l  _  0  0  0  4”  -42)  421  1 1, 1  (27) 

P2  ~  i<»  C»>  d{"  0  0  0  1  P2 

t  *P  4"  4°  0  0  0  Z. 

U"  4"  4"  0  0  0  j 
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For  convenience  in  the  following  derivations,  let  us  define 


T!_2  =  I  tl 

1  P2 


P? 

s? 

•  T°°  = 

’  ’-2  _  I  p r 

<t°° 

s2 

t°° 

l2 


;  t,_2  = 


and  cx  =  the  6x6  matrix  in  Eq.  (27).  With  these  notations  at  hand,  Eq.  (27)  can  be 
rewritten  as 

Ti_2  =  0:-Ti_2  (29) 

Since  Tj_2  =  T“  2  +  Ti_2,  Tj_2  can  be  solved  from  Eq.  (29) 

Ti_2  =  K  •  T“  2  !  where  K  =  a  •  (I  —  a)-1  (30) 

Further,  we  define  T  =  (pi,  Si,  t\)T  and  Kj  =  the  first  three  rows  of  K.  Hence,  from  Eq.  (30), 
the  perturbed  stresses  on  surfaces  of  a  microcrack  due  to  the  existence  of  a  second  microcrack 
are: 

T  =  K1-T“2  (31) 

Finally,  since  the  two  microcracks  under  consideration  are  aligned,  it  can  be  easily  shown 


rpoo  _ 

A  1—2  — 


p?) 

0 

0 

1 

0 

0 

sl° 

0 

0 

0 

0 

0 

tr  _ 

0 

0 

0 

0 

1 

P2  ~ 

0 

0 

1 

0 

0 

c°° 

*2 

0 

0 

0 

0 

0 

400 

l2 

0 

0 

0 

0 

1 

=  K,  .  T° 


From  Eq.  (32)  and  (31),  we  conclude  that  the  average  perturbed  stress  T  over  surfaces  of  a 
microcrack  is  simply: 

T  =  K  •  t°°  ;  where  KeK,  K2  (33) 

Substitution  of  Eq.  (33)  into  (13)  then  renders  realizations  for  (T)  in  the  previous  section. 
Therefore,  the  ensemble  average  approach  is  completely  defined.  For  non-aligned  penny¬ 
shaped  microcracks,  more  complicated  (though  analogous)  derivations  will  be  involved. 


III.2.3.  Some  test  problems  for  two-microcrack  interaction 


A  number  of  test  problems  are  considered  in  this  section  to  examine  the  performance 
of  the  approximate  analytical  solutions  presented  in  Sec.  2.2.  These  include  normal  and 
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shear  loadings  for  two  aligned  coplanar  or  stacked  penny-shaped  microcracks.  It  is  not 
our  intention,  however,  to  propose  highly  accurate  analytical  solutions  to  compute  local 
stresses  and  stress  intensity  factors  (SIFs)  at  all  points  on  microcrack  surfaces  for  a  two- 
microcrack  interaction  problem.  Instead,  reasonably  accurate  analytical  average  perturbed 
stresses  (in  terms  of  elementary  functions)  over  microcrack  surfaces  are  sought  in  order  to 
exploit  the  ensemble  average  approach.  In  fact,  if  one  is  interested  in  deterministic  3-D 
microcrack  interaction,  excellent  numerical  method  has  been  proposed  by  Kachanov  and 
Laures  (!989).  It  is  noted  that  key  steps  in  Kachanov  and  Laures  (1989)  method  also  focus 
on  the  computation  of  average  pairwise  “transmission  factors”  and  average  tractions.  Once 
average  tractions  become  known,  one  can  certainly  compute  projected  local  stresses  and  SIFs 
at  any  point  on  microcrack  surfaces.  In  general,  the  simple  analytical  solutions  presented  in 
Sec.  2.2  are  not  as  accurate  as  those  proposed  in  Kachanov  and  Laures  (1989).  Nonetheless, 
the  latter  relies  on  extensive  numerical  computations  of  transmission  factors  for  all  points 
on  microcrack  surfaces  and  is  therefore  not  employed  here. 

It  is  emphasized  that,  given  an  existing  microcrack,  the  location  of  the  second  microcrack 
is  random  within  the  ensemble  average  framework.  Although  stress  interactions  between 
two  moderately  spaced  random  microcracks  are  not  very  strong,  their  cumulative  effects  are 
important  due  to  high  spatial  probability.  On  the  other  hand,  stress  interactions  between  two 
closely  spaced  random  microcracks  are  strong  yet  their  spatial  probability  is  lower.  Therefore, 
contributions  from  both  closely  and  moderately  spaced  microcracks  should  be  accounted  for 
in  the  ensemble  average  approach.  The  ensemble-average  spacing  of  microcrack  arrays, 
clearly,  depends  on  microcrack  concentration. 

Case  I:  Two  equal-size  coplanar  microcracks  under  normal  loading.  Though  numerical 
results  for  SIFs  are  available  for  two  coplanar  microcracks  under  normal  and  shear  loadings 
(Fabrikant  (1987,  1989)),  exact  results  for  average  tractions  projected  over  microcrack  sur¬ 
faces  are  not  documented.  Nevertheless,  Kachanov  and  Laures  (1989)  show  that  their  results 
for  SIFs  are  very  close  to  those  of  Fabrikant  (1987,  1989)  for  two  equal-size  coplanar  micro¬ 
cracks.  Therefore,  the  average  tractions  computed  by  Kachanov  and  Laures  (1989)  should 
be  quite  accurate  in  the  case  of  coplanar  microcracks.  The  ratios  of  perturbed  vs.  far-field 
average  normal  stresses  (p/p00)  are  listed  in  Table  1  for  various  l /2a  values  (the  Poisson 
ratio  v  —  0.25).  Here,  l  signifies  the  center-to-center  distance  between  two  microcracks  and 
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2 a  is  the  microcrack  size.  For  example,  l /2a  =  1.00025  means  that  the  smallest  distance  be¬ 
tween  microcrack  tips  is  only  0.00025a.  The  results  reported  in  Kachanov  and  Laures  (1989) 
and  relative  differences  between  the  two  results  are  also  given  in  Table  1  for  comparison. 
From  Table  1,  it  is  clear  that  interaction  renders  average  stress  amplification.  The  effect  of 
microcrack  interaction  decays  as  the  distance  between  two  microcracks  increases.  Moreover, 
the  differences  between  our  simple  calculations  and  those  of  Kachanov  and  Laures  (1989) 
are  very  small  in  general. 

Case  II:  Two  equal-size  coplanar  microcracks  under  shear  loading.  The  ratios  of  per¬ 
turbed  vs.  far-field  average  shear  stresses  (t/t00)  are  listed  in  Table  2  for  various  //2a 
values  (the  Poisson  ratio  v  —  0.5).  From  Table  2,  it  is  seen  again  that  interaction  renders 
average  stress  amplification.  The  differences  between  our  simple  calculations  and  those  of 
Kachanov  and  Laures  (1989)  are  again  rather  small;  see  Table  2. 

Case  III:  Two  equal-size  stacked  microcracks  under  normal  loading.  The  ratios  of  per¬ 
turbed  vs.  far-field  average  normal  stresses  (p/p°°)  are  listed  in  Table  3  for  various  //2a 
values  (the  Poisson  ratio  v  =  0.25).  Obviously,  from  Table  3,  interaction  renders  strong 
average  stress  shielding.  The  degree  of  shielding  in  this  case  is  much  stronger  than  the  degree 
of  amplification  in  Cases  I  or  II.  Fabrikant  (1989)  did  not  provide  numerical  results  for  SIFs 
nor  average  stresses  for  two  stacked  microcracks  under  normal  or  shear  loadings.  Therefore, 
exact  results  are  not  available.  For  comparison,  the  differences  between  our  simple  calcula¬ 
tions  and  those  of  Kachanov  and  Laures  (1989)  are  listed  in  Table  3.  It  is  observed  that  the 
difference  is  not  very  significant  if  //2a  is  greater  than  0.5.  The  difference  increases  as  the  two 
microcracks  move  closer,  due  to  sharp  variations  of  stress  fields  in  the  close  neighborhood  of 
inter*  'ting  microcracks. 

Case  IV:  Two  equal-size  stacked  microcracks  under  shear  loading.  The  ratios  of  per¬ 
turbed  vs.  far-field  average  normal  stresses  (p/p°°)  are  given  in  Table  4  for  various  //2a 
values  (the  Poisson  ratio  v  =  0.25).  From  Table  4,  interaction  renders  average  stress  shield¬ 
ing  when  the  distance  /  is  small  and  very  weak  stress  amplification  when  //2a  >  0.5.  The 
differences  between  our  simple  calculations  and  those  of  Kachanov  and  Laures  (1989)  are 
small  for  //2a  greater  than  0.25;  see  Table  4. 

Since  pairwise  microcrack  interactions  are  random  (probabilistic)  within  the  framework 
of  the  ensemble  average  approach,  the  errors  associated  with  the  present  approximate  anal¬ 
ysis  should  be  statistically  averaged  over  all  possible  realizations.  Therefore,  the  (pairwise) 
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ensemble-average  error  of  the  present  method  should  be  small  as  long  as  the  microcrack 
concentration  is  not  too  high. 
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With  the  ensemble  average  framework  and  analytical  pairwise  interaction  solutions  at 
hand,  we  are  now  ready  to  construct  the  ensemble- average  constitutive  equations  of  brittle 
solids  with  many  randomly  located,  interacting  microcracks.  By  substituting  Eq.  (10)  and 
(33)  into  (9)  and  taking  the  ensemble  average,  we  arrive  at 


(T)  =  (K0  +  /(x)(K»  ■  r“ 


(34) 


where 


K  dX! 


J  Kr2  sin  rp  drdtpdO 


(35) 


In  the  foregoing  equations,  it  is  implicitly  assumed  that  all  microcracks  are  of  equal  size  and 
of  same  orientation.  In  addition,  the  spherical  coordinate  system  (r,ip,0)  is  used  to  describe 
the  random  location  (xj)  of  the  second  microcrack  relative  to  the  first  random  microcrack 
centered  at  x.  Note  that  i/>  varies  from  0  to  it  and  0  ranges  from  0  to  2ir.  If  we  normalize  r 
with  respect  to  the  microcrack  radius  a  (i.e.,  £  =  r/a),  then  Eq.  (35)  can  be  recast  as 


(K)=<z3(K)  =  <.3jf 


Kf  2  sin  rj>  d£dipd0 


(36) 


Combining  Eq.  (6),  (34)  and  (36),  we  obtain  the  local  ensemble- average  damage-induced 
strain  (at  a  typical  point  x): 

(0(x)  =  (S*)(x)  •  t°°  (37) 

Here  the  ensemble  averaged,  damage-induced  local  compliance  has  two  components: 

(S')(x)  =  (S*1)(x)  +  (S‘2)(x)  (38) 


where 

(S-’)(x)  =  ~  ~r /(x)a3g  ■  K„ 

3L(2  ~  '  (39) 

<S'I)(X)  =  ^2T7)/!(x)°6g  • 

It  is  noteworthy  that  Eq.  (39)  actually  reveals  physical  nonlocal  effects  in  constitutive  equa¬ 
tions  through  the  ensemble  averaging  process.  That  is,  the  stress-strain  laws  at  a  material 
point  x  within  a  RVE  depend  on  the  constitutive  laws  of  all  neighboring  points.  This  is  a 
physical  nonlocal  approach,  at  variance  with  postulated  nonlocal  theories  due  to  Eringen  and 
Edelen  (1972). 
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To  obtain  volume  averaged  moduli  due  to  microcracks  within  a  RVE,  one  simply  applies 
the  volume-average  operator  to  Eq.  (38)-(39).  As  a  consequence,  we  have 

W)  =  i  /(S-X*>  *  =  ^  [jf  <S’')(x)  +  ^(Srt)(x)  <6c]  (40) 


where 


ToTiT  16(1  —  I/2)  v  fvf(x)dx 
{s  >=3£(2-,)g-K““  - V - 

(S^y  =  HI  .zn’jg .  (K) 

'  '  3£(2-«/)g  '  '  V 

Let  us  assume  that  there  are  N  microcracks  in  the  RVE;  i.e, 


f(x)  dx  = 


Further,  consider  the  case  in  which  the  variance  of  the  PDF  (/(x))  for  locations  of  microc¬ 
racks  is  small  (e.g.,  uniform  probability).  We  may  therefore  write 


[/v,  /(x)  dx]  JV2 


Substitution  of  Eq.  (42)  and  (43)  intc  (41)  then  leads  to 

7  ■  -r  16(1  —  i/2) 
j—  16(1 -i/2)  -  2 

<S  >=3£(2-^)g'(K>W 

where 


is  the  (volume-average)  microcrack  concentration  parameter.  It  is  noted  that  (S*1)  actually 
corresponds  to  the  first  order  contribution  due  to  non-interacting  microcracks;  i.e.,  the  simple 
Taylor’s  model  is  recovered.  Moreover,  (S*2)  represents  the  second  order  contribution  due 
to  pairwise  microcrack  interaction. 

Finally,  the  overall  (volume-ensemble  averaged)  effective  moduli  for  a  microcrack  weak¬ 
ened  solid  is  obtained  by  adding  the  elastic  compliance  S°  to  (S): 


(S)  =  S°+(S*1)  +  (S*2) 


(46) 
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In  particular,  in  the  case  of  3-D  linear  isotropic  elasticity,  S°  reads 

1  — v  — v  0  0  0 

-v  1  -v  0  0  0 

-v  —v  1  0  0  0 

0  0  0  2(1  +  //)  0  0 

0  0  0  0  2(1  +  //)  0 
0  0  0  0  0  2(1  +  //) 

It  is  emphasized  that  the  local  macroscopic  constitutive  laws  for  a  RVE  are  recovered  by  the 
volume  averaging  process. 
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In  this  section,  some  numerical  examples  involving  aligned  microcracks  are  presented  to 
illustrate  the  proposed  ensemble- volume  average  approach  in  three-dimension.  In  the  first 
example,  effects  of  microcrack  interactions  are  neglected  and  the  well-known  Taylor’s  model 
is  recovered.  Subsequently,  the  proposed  second  order  (in  microcrack  concentration)  micro- 
crack  interaction  model  is  implemented.  Finally,  we  compare  the  results  of  the  present  ap¬ 
proach  with  some  existing  methods,  including  the  Taylor’s  model,  the  self-consistent  method 
and  the  differential  scheme.  Closed  microcrack  contributions  to  compliances  are  neglected 
in  this  work. 


III.4.1.  Dilute  non-interacting  aligned  microcracks 

Let  us  consider  the  case  in  which  effects  of  microcrack  interactions  are  totally  neglected. 
This  can  be  done  by  simply  dropping  the  term  (S’2)  in  Eq.  (46): 

Jsj  =  s°  +  WT)  (48) 


From  Eq.  (8),  (11),  and  (*4),  we  obtain 


(S*1)  = 


'0  0  0  0  0  0 

„  0  0  0  0  0  0 

16(1  -  v2)u  0  0  2  -  1/  0  0  0 

3£(2  -i/)  0  0  0  000 

0  0  0  0  2  0 

0  0  0  0  0  2 


Combination  of  Eq.  (47)  and  (49)  then  yields 


'  CO  CO  CO 

On  Oj2  J13 

50  CO  CO 

21  ^22 


<S)  = 


SO  C  0 

3i  J3: 


’32  ‘-’33 

0  0 


o  s°4 


where  1  =  xx,  2  =  yy,  3  =  zz,  4  =  xy,  5  =  yz,  6  =  zx  (Voigt’s  notation);  5®  components 
are  given  in  (47),  and 

533  =  5^  +  5^ 

sis  =  s°K  +  s;l  (5i) 

See  =  Sj,  + 
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in  which 


.1  _  16(1  -  v2) 


c*i  = 
‘-’33  — 


‘-’56 


3  E 


u 


_  c.1  =  32(1  —  v2) 
66  “  3£(2  -  v) 


u> 


(52) 


It  is  observed  that  only  the  three  compliance  components  in  Eq.  (51)  are  changed  due 
to  the  presence  of  randomly  located,  aligned,  penny-shaped  microcracks.  One  notes  that 
Eq.  (50)-(52)  indeed  recover  the  well-known  Taylor’s  model.  Further,  the  effective  constitu¬ 
tive  equations  belong  to  the  category  of  transverse  isotropy.  This  is  typical  of  fiber  breaks  in 
unidirectional  reinforced  fiber  composites. 


III.4.2.  Aligned  interacting  penny-shaped  microcracks 

We  now  focus  on  the  proposed  ensemble- volume  averaged,  pairwise  interaction  model. 
For  a  given  microcrack  concentration  u>,  the  overall  effective  moduli  in  Eq.  (46)  can  be  eval¬ 
uated  by  carrying  out  the  integration  in  Eq.  (36).  That  is,  one  needs  to  perform  integration 
for  the  matrix  K  over  the  active  microcrack  domain  (if,  ip,  6).  It  is  emphasized  that  the 
values  of  K  decay  rapidly  as  the  distance  between  a  pair  of  random  microcracks  increase. 
Therefore,  compliance  contributions  due  to  remote  integration  region  can  be  neglected. 

The  integration  in  Eq.  (DC)  can  be  effectively  computed  by  the  Gauss  quadrature  scheme 
with  three  independent  variables  -  £,  ip  and  0;  see  Table  5  for  convergence  behavior.  Our 
numerical  experiments  show  that  use  of  fmax  =  20  (or  rmax  =  20a)  is  quite  acceptable;  i.e., 
contributions  from  the  £  >  20  domain  can  indeed  be  neglected.  In  Table  5,  the  minimum 
radius  of  integration  for  r  (or  f)  is  assumed  to  be  rmi„  =  2a  (or  £min  =  2).  The  only 
nonzero  components  in  (K)  are  simply  (X'13),  (A'26)  and  (A35)  (the  latter  two  are  equal).  It 
is  observed  that  the  normal  component  ((A'13))  is  much  greater  than  the  shear  components 
((A^)  and  (A'35).  Furthermore,  use  of  rmin  =  2a  implies  that  the  minimum  allowable 
distance  between  any  two  microcrack  centers  is  2a.  The  '  fore,  according  to  the  “face-center 
cubic”  calculation,  the  maximum  allowable  u  is  1/4 \/2  ~  0.1768  for  rmin  =  2a.  Note  that 
the  maximum  allowable  microcrack  density  would  be  0.741  if  u>  is  defined  as  4nNa3/3V 
(Budiansky  and  O’Connell,  1976). 

To  achieve  higher  maximum  permissible  microcrack  densities,  various  minimum  radii  of 
integration  rmtn  are  used  to  compute  (K);  see  Table  6  for  details.  The  value  of  fmax  is 
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chosen  as  20.0  in  Table  6.  We  observe  significant  increases  in  both  (K13)  and  (A'26)  as 
£min  decreases.  The  maximum  microcrack  concentration  u>  corresponding  to  the  £mtn  =  1.1a 
case  is  found  to  be  1.0627  (or  4.451  if  the  definition  in  Budiansky  and  O’Connel]  (1976)  is 
employed).  It  is  emphasized  that  Table  6  is  not  a  summary  of  convergence  behavior,  but  a 
display  of  (K)  for  different  maximum  allowable  microcrack  densities. 

A  A  A  A  A 

For  convenience,  let  us  define  (K 13)  =  ki  and  (K2 e)  =  (A35)  =  k2.  Therefore,  according 
to  Eq.  (8)  and  (44),  we  obtain 


16(1  -  *»)  a 
3£(2  -  v) 


'0  0  0  0 

0  0  0  0 

0  0  (2-v)kx  0 

0  0  0  0 

0  0  0  0 

.0  0  0  0 


0 

0 

0 

0 

2  k2 
0 


(53) 


Adding  Eq.  (53)  to  (50),  we  thus  arrive  at  the  expression  for  overall  effective  complaince  with 
second  order  microcrack  interaction.  Apart  form  the  three  compliance  components  (S33), 
( S55 )  and  (S^),  it  is  seen  that  all  other  components  are  identical  to  the  elastic  components. 
In  particular,  the  three  compliance  components  which  change  under  microcrack  interaction 
take  the  following  explicit  forms: 


(S33)  =  S%3  H - — -  (w  +  kjut2) 

(S55)  =  (See)  =  S°5  +  3^2 '-V)  ^  +  k2u2) 


III. 4. 3.  Comparison  with  some  existing  methods 

The  effects  of  fiber  breaks  and  aligned  penny-shaped  microcracks  on  the  stiffness  of 
unidirectional  fiber  composites  have  been  studied  extensively  in  the  literature.  We  refer 
to  Laws  and  Dvorak  (1987)  for  excellent  presentation  by  using  the  self-consistent  method 
and  the  differential  scheme.  It  is  noted  that,  in  Laws  and  Dvorak  (1987),  the  microcrack 
density  is  defined  as  a  =  8Na3/V ;  i.e.,  their  o  is  equal  to  8u  in  this  paper.  The  maximum 
microcrack  concentration  considered  in  Laws  and  Dvorak  (1987)  is  u>  =  a/8  =  0.125.  The 
macroscopic  (overall)  material  behavior  is,  not  surprisingly,  transversely  isotropic.  It  should 
be  realized  that  the  Taylor’s  model,  the  self-consistent  method  and  the  differential  scheme 
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all  belong  to  the  category  of  “effective  medium”  approaches.  Namely,  theses  methods  do  not 
depend  on  the  distributions  of  microcrack  locations  at  all. 

In  the  following  numerical  computations,  the  Young’s  modulus  E0  is  taken  as  0.5  MPa 
and  the  shear  modulus  Go  as  0.2  MPa  ( v  =  0.25)  for  the  virgin  matrix  material.  The  overall 
(ensemble- volume  averaged)  longitudinal  normal  and  shear  compliances,  (S33)  and  (S55),  are 
plotted  against  the  microcrack  concentration  parameter  a?  in  Figures  2  and  3,  respectively. 
The  values  of  k\  and  k 2  in  Eq.  (54)  are  taken  as  1.049792  and  0.355988,  respectively,  which 
correspond  to  =  1.1  (i.e.,  the  maximum  allowable  density  u  =  1.0627).  Note  that,  once 
more,  u  —  1  corresponds  to  a  =  8  in  Laws  and  Dvorak  (1987).  For  comparison,  the  results 
obtained  by  using  the  Taylor’s  model,  the  self-consistent  method  and  the  differential  scheme 
are  also  displayed  in  Figures  2  and  3.  We  would  like  to  comment  that:  (a)  the  Taylor’s 
model  is  really  suitable  for  dilute  microcrack  concentrations;  (b)  the  self-consistent  method 
and  the  differential  scheme  are  suitable  for  low  or  moderate  u>\  and  (c)  the  proposed  statistical 
pairwise  microcrack  interaction  model  is  suitable  for  moderately  high  (not  extremely  high) 

U). 


Figures  4  and  5  show  the  effect  of  microcrack  density  on  the  normalized  longitudinal 
Young’s  modulus  El/ Eq  and  shear  modulus  Gl/Gq  for  four  different  models.  We  have  em¬ 
ployed  the  standard  notation:  El  =  ( S33 )  and  Gl  =  (S55).  In  addition,  it  is  straightforward 
to  express  the  normalized  moduli  in  Figures  4  and  5  as  follows  for  the  Taylor’s  model 


El  1 

Eq  1  -}-  5  U! 


Gl  1 

G0  1  +  2.28571u; 


and  for  the  present  model 


El  1 

Eq  1  +  5u;(l  +  1.05w) 


Gl  = _ 1 _ 

G0  1  -(-  2.2857lu>(l  +  0.356u>) 


(55) 


(56) 
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III.5.  Higher  order  ensemble-average  formulation  of  microcrack  in¬ 
teraction 


The  statistical  micromechanical  model  presented  so  far  is  based  on  the  concept  of  pair¬ 
wise  microcrack  interaction.  This  pairwise  microcrack  interaction  mechanism  essentially 
corresponds  to  a  second  order  damage  theory  since  it  involves  terms  of  u>2.  The  classical 
Taylor’s  model,  on  the  other  hand,  is  a  first  order  microcrack  theory  since  it  involves  only 
terms  of  u.  Within  the  context  of  the  ensemble-volume  average  approach,  one  can  system¬ 
atically  incorporate  many-microcrack  interaction  mechanisms  into  the  proposed  ensemble 
average  framework  (see  Ju  and  Chen  (1990.  .  In  essence,  one  needs  to  re-derive  the  en¬ 

semble  average  of  the  perturbation  in  local  stress  field  due  to  n-microcrack  interactions 
(n  >  3);  see  Eq.  (9)-(12).  That  is,  by  formulating  and  computing  higher  order  corrections  on 
local  stress  fields  (T)  due  to  neighboring  microcrack  interactions,  one  can  construct  higher 
oder  microcrack  interaction  models. 


To  illustrate  the  foregoing  statements,  let  us  consider  a  three-microcrack  (third  order  in 
u)  interaction  mechanism  within  the  proposed  framework  of  the  ensemble-volume  average 
approach.  For  simplicity,  we  assume  that  all  microcracks  are  aligned  and  of  equal  size. 
Following  the  definitions  and  assumptions  described  in  Section  2.1,  one  can  lecapitulate  the 
local  ensemble  stress  perturbation  in  Eq.  (9): 

<T) +  <*>  =  <{  |  }>  +  < 

where  (T)  is  the  first  order  local  ensemble  stress  perturbation  due  to  pairwise  microcrack 
interactions  (see  Eq.  (9),  (12)),  and  (T)  is  the  second  order  local  ensemble  stress  perturbation 
due  to  the  higher  (third)  order  microcrack  interactions.  In  particular,  (T)  can  be  expressed 
as  (cf.  Eq.  (12)) 

(i’)  =  ^{T)(x;xj|x2)/(x2|x;x,)  dx2  (58) 

Here,  (T)(x;xijx2)  is  the  second  order  ensemble- average  stress  perturbation  of  a  microcrack 
centered  at  x,  given  a  microcrack  centered  at  Xi,  over  a  subclass  of  reali  sations  which  have  a 
microcrack  centered  at  X2.  Further,  /(x2(x;  Xi)  is  the  conditional  probability  density  function 
(PDF)  for  finding  a  microcrack  centered  at  x2  given  two  microcracks  fixed  at  x  and  at  Xj. 
The  conditional  PDF  /(x2|x;Xj)  can  be  further  simplified  to  /(x)  by  the  assumptions  of  local 


|J)  (57) 
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homogeneity  and  reasonable  randomness  (i.e.  statistical  independence).  The  active  (open) 
integration  domain  s  depends  on  loading  conditions  and  ranges  of  microcrack  interactions. 

The  solutions  of  T  (or  the  total  stress  vector  T)  for  a  system  of  three  (or  many)  arbi¬ 
trarily  located  microcracks  were  previously  investigated  by  Kachanov  and  Laures  (1989).  In 
particular,  Kachanov  and  Laures  (1989)  pursued  extensive  numerical  computations  to  ob¬ 
tain  the  “transmission  factors”  for  local  stresses.  Alternatively,  one  can  derive  approximate 
analytical  closed-form  solutions  for  T  by  following  the  procedures  presented  in  Sec.  2.2  of 
this  paper. 

For  clarity,  let  us  express  T  and  (e*)(x)  as  follows  (cf.  Eq.  (34),  (37)  and  (38),  assuming 
local  homogeneity  and  reasonable  randomness): 


> 

(T)  =  (T~  +  T  +  T)  =  (K0  +  /(x)(K)  +  /2(x)(K'))  •  r°° 

(59) 

(e*)(x)  =  {(S*1)(x)  +  (S*2)(x)  +  (S*3)(x)}  •  r°°  =  (S*)(x)  ■  r°° 

(60) 

where  T  =  K'  r°°  and  (cf.  Eq.  (39)) 

i 

(61) 

(62) 

> 

<s'3><x>=^,Wv<^ 

(63) 

In  Eq.  (63),  we  have  defined  (K')  =  1/a6  (K').  It  is  emphasized  that  both  K0 

and  K  are 

expressed  explicitly  in  closed-form  formulas.  Similarly,  K'  (or  t)  can  also  be  constructed 
in  closed-form  as  follows.  One  starts  by  expanding  Eq.  (14)  into  nine  linear  equations  with 
j  =  1,2,3.  Then,  one  obtains  expressions  similar  to  Eq.  (22),  with  the  understanding  that 
permutations  1-2,  2-3,  3-1  are  involved.  Eq.  (24)  in  Sec.  2.2  is  modified  to  include  the  third 
microcrack’s  contribution  to  stress  perturbations.  Subsequently,  Eq.  (27)  is  expanded  to  a 
9  by  9  system  with  a  denoting  a  9  by  9  coefficient  matrix.  Therefore,  we  arrive  at  explicit 
formulas  similar  to  Eq.  (30)-(33);  and  K',  T  can  be  expressed  in  closed-form  as  in  the  two- 
microcrack  interaction  problem.  Finally,  it  is  noted  that  the  computation  of  (K')  involves 
integration  over  the  domain  of  all  possible  positions  of  two  active  neighboring  microcracks 
(cf.  Eq.  (36)). 
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Following  the  standard  procedures  presented  in  Sec.  3,  it  can  be  shown  that  (S*3)  in 
Eq.  (39)  introduces  the  third  order  terms  (in  u>3)  to  overall  compliance  due  to  third  order 
microcrack  interactions.  Therefore,  a  third  order  statistical  micromechanical  model  can  be 
actually  constructed.  By  repeating  the  foregoing  procedures,  we  can  formulate  a  complete 
(though  very  complicated)  hierarchical  family  of  statistical  microcrack  theories  of  arbitrary 
(desired)  orders. 
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III.6.  Conclusions 

An  innovative  three-dimensional  statistical  micromechanical  theory  for  microcrack-weakened 
brittle  solids  is  presented  based  on  the  concepts  of  ensemble-volume  average  and  pairwise 
microcrack  interaction.  A  physically  “ nonlocal ”  description  of  the  material  behavior  can 
be  obtained  during  the  ensemble  average  process.  The  overall  compliances  of  microcrack- 
weakened  brittle  solids  are  derived  by  performing  ensemble- volume  averaged  integration  over 
the  entire  domain  of  a  representative  unit  cell.  To  account  for  microcrack  interaction  effects, 
the  ensemble  average  integration  is  performed  over  a  finite  “interaction  radius”.  The  pro¬ 
posed  approach  is  fundamentally  different  from  existing  effective  medium  methods  which  do 
not  depend  on  locations  and  configurations  of  microcracks,  Further,  the  proposed  microcrack 
interaction  framework  does  not  require  the  use  of  Monte  Carlo  simulation.  Some  numer  cal 
examples  are  given  to  illustrate  the  proposed  model.  The  resulting  predictions  are  compared 
with  other  existing  methods.  Higher  order  microcrack  interaction  formulation  is  also  briefly 
presented. 
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Table  1.  Two  equal-size  coplanar  microcracks  under  normal  loading  {v  =  0.25) 


//2a 

Present 

Kachanov 

Difference(%) 

1.00025 

1.0354 

1.0837 

4.46 

1.005 

1.0348 

1.0779 

4.00 

1.05 

1.0296 

1.0529 

2.21 

1.1 

1.0251 

1.0398 

1.41 

1.15 

1.0215 

1.0315 

0.97 

1.25 

1.0161 

1.0214 

0.52 

1.5 

1.0088 

1.0104 

0.16 

2.0 

1.0035 

1.0038 

0.03 

2.5 

1.0018 

1.0019 

0.01 
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Table  2.  Two  equal-size  coplanar  microcracks  under  shear  loading  {v  =  0.5) 


J/2a 

Present 

Kachanov 

Difference(%) 

1.005 

1.0681 

1.1017 

3.05 

1.05 

1.0580 

1.0703 

1.15 

1.25 

1.0317 

1.0292 

0.24 

1.5 

1.0174 

1.0144 

0.30 

1.75 

1.0106 

1.0084 

0.22 

2.0 

1.0070 

1.0054 

0.16 

2.5 

1.0035 

1.0026 

0.09 

3.5 

1.0013 

1.0006 

0.07 
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Table  3.  Two  equal-size  stacked  microcracks  under  normal  loading  ( v  =  0.25) 


//2a 

Present 

Kachanov 

Difference(%) 

0.05 

0.5004 

0.5583 

10.37 

0.25 

0.5383 

0.6689 

19.52 

0.35 

0.5836 

0.7158 

18.47 

0.5 

0.6667 

0.7777 

14.27 

0.75 

0.7928 

0.8562 

7.40 

1.0 

0.8754 

0.9073 

3.52 

1.5 

0.9505 

0.9588 

0.87 
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Table  4.  Two  equal-size  stacked  microcracks  under  shear  loading  ( v  =  0.25) 


//2a 

Present 

Kachanov 

Difference(%) 

0.05 

0.5549 

0.6613 

16.09 

0.25 

0.8214 

0.8886 

7.56 

0.5 

1.0002 

0.9837 

1.68 

0.75 

1.0233 

1.0053 

1.79 

1.0 

1.0180 

1.0084 

0.95 

1.25 

1.0120 

1.0072 

0.48 

1.5 

1.0080 

1.0055 

0.25 

2.0 

1.0039 

1.0031 

0.08 

2.5 

1.0021 

1.0026 

0.05 
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Table  5.  Convergence  behavior  of  (K) 


£max  —  (-^13) 


4.0 

0.150807 

6.0 

0.168554 

8.0 

0.173029 

10.0 

0.174644 

20.0 

0.176139 

40.0 

0.176327 

80.0 

0.176351 

160.0 

0.176354 

320.0 

0.176354 

640.0 

0.176354 

VS.  tmax  (given  £mtn  =  2) 

(*»>  =  (tf35> 

0.059503 

0.066656 

0.068467 

0.069121 

0.069728 

0.069804 

0.069828 

0.069829 

0.069829 

0.069829 
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Table  6.  Numerical  integration  of  (K)  for  different  £mtn  (given 


.  —  f'mt'n/® 

(Kts) 

II 

c* 

2.0 

0.176139 

0.069728 

1.75 

0.305726 

0.105974 

1.5 

0.601154 

0.195689 

1.25 

0.850189 

0.281202 

1.125 

0.975302 

0.327491 

1.1 

1.049792 

0.355988 

=  20) 
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The  parameters  6,,  c,-  and  di  in  Eq.  (22)  can  be  shown  to  be: 


,  2  .  az 2 

b\  =  —  (1  +  2v)k\  -i - k2 

*  92 

2  J  z2 

hi  =  —ag£  (1  —  2v)k\  -f  — k3  cos  2$ 

*  92  J 

2  az 2 

b 3  =  -  *X - fc2 

*  92  \ 

2  ,  z2 

fe4  =  — a<73  (1  —  2i/)A:j  -( - k3  sin  20 

v  9  2  . 

65  - - zg29zkA  sin  0 

7T 

b6  =  --zg2g3k4cos<f> 

7T 


d  =  /3  COS  0 

c2  =  /1  cos  0  +  /5  cos  30  +  /6(cos  <j>  +  cos  30) 
c3  =  /4  sin  <f>  +  /5  sin  30  +  /6(sin  30  +  sin  0) 
c4  =  — 2/gCos  0 
c5  =  /7  +  /s  COS  20 
ce  =  fs  sin  20 


di  =  fz  sin  0 

d2  =  —fA  sin  <f>  +  fs  sin  30  +  /6(sin  30  —  sin  0) 
d3  =  fA  cos  0  —  /5  cos  30  +  fe{cos  0  —  cos  30) 
dA  =  — 2/6sin  0 
ds  =  /s  sin  20 


d6  =  fr  -  fs  cos  20 


(66) 
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where 


/i  = 
/2  = 


<h 

9s 

92 

9s 


h  = 


8 


?r(2  —  V ) 


[-2(1  +  f)ag3f,  +  zgigikt] 


,  8(1--).., 
u  =  -7^7)aS3!' 

h  =  -W=T)Z93hs< 
h  =  7WlT)Z9’93k' 

!l  =  l  (*'  + 


fa  = 


tt(2  —  v ) 


sK-a/i  +  zjife) 


=  („’  - /2)>'2 
,2  =  (/2  -  a2)1'1 

,34 

h 

a 

9*  =  T 

9s  =  ll-ll 

*i  =  a/2  -  sin-1  (#4) 

.  l\  +  a2(2a2  +  2z2  -  3 p2) 

k2  =  - 3 - 

9s 

,  a2(6/2  -  21]  +  p 2)  -  5/3 

fc3  =  - 3 - 

9s 

,  a2(4 1\  -  5 ,2)  +  /< 

*4  = - -3 - 


(67) 


(68) 


(69) 
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Figure  1. 
Figure  2. 
Figure  3. 
Figure  4. 
Figure  5. 


Coordinate  systems  for  two  microcracks. 

Comparison  of  overall  longitudinal  normal  compliance  (533). 
Comparison  of  overall  longitudinal  shear  compliance  (S55). 
Comparison  of  normalized  longitudinal  Young’s  modulus. 
Comparison  of  normalized  longitudinal  shear  modulus. 


(S33l  (MPa) 


(Sss)  (MPa) 


Normalized  Young’s  modulus  EL/E0 


Microcrack  concentration  u 


Figure  4 


Comparison  of  normalized  longitudinal  Young’s  modulus. 


Normalized  shear  modulus  Gl/G 


